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Abstract. Experimental Physics with Relativistic Heavy Ions dates from 1992 when 
a beam of ^^^Au of energy greater than lOA GeV/c first became avalilable at the 
Alternating Gradient Synchrotron (AGS) at Brookhaven National Laboratory (BNL) 
soon followed in 1994 by a ^°^Pb beam of 158A GeV/c at the Super Proton Synchrotron 
(SPS) at CERN (European Center for Nuclear Research). Previous pioneering 
measurements at the Berkeley Bevalac ^ in the late 1970's and early 1980's were 
at much lower bombarding energies (~ lA GeV/c) where nuclear breakup rather than 
particle production is the dominant inelastic process in A+A collisions. More recently, 
starting in 2000, the Relativistic Heavy Ion Collider (RHIC) at BNL has produced 
head-on collisions of two lOOA GeV beams of fully stripped An ions, corresponding 
to nucleon-nucleon center-of-mass energy, ^/sWn = 200 GeV, total cm. energy 200A 
GeV. The objective of this research program is to produce nuclear matter with extreme 
density and temperature, possibly resulting in a state of matter where the quarks and 
gluons normally confined inside individual nucleons (r < Ifm) are free to act over 
distances an order of magnitude larger. Progress from the period 1992 to the present 
will be reviewed, with reference to previous results from light ion and proton-proton 
collisions where appropriate. Emphasis will be placed on the measurements which 
formed the basis for the announcements by the two major laboratories: "A new state 
of matter" , by CERN on Feb 10, 2000 and "The perfect fluid" , by BNL on April 19, 
2005. 
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1. Introduction 

In the early 1970's, it became clear that the nucleon was not an elementary particle 
but was composed of a substructure of 3 valence quarks confined into a bound state 
by a strong interaction, Quantum Chromo Dynamics (QCD), which is mediated by the 
exchange of color-charged vector gluons. j2] In sharp distinction to the behavior of the 
uncharged quanta of quantum electrodynamics (QED), the color-charged gluons of QCD 
interact with each other. This leads to the property of asymptotic freedom, [3 |3] the 
reduction of the effective coupling constant at short distances, and is believed to provide 
the confinement property at long distances where the quarks and gluons behave as if 
attached to each other by a color string. It is worth reviewing some developments leading 
up to and immediately following the discovery of QCD which bear on the importance 
of Relativistic Heavy Ion collisions as a probe of this fundamental interaction. 

In the early 1960 's, with the construction of proton accelerators with energies 
well above the threshold for anti-proton production, a veritable 'zoo' of new particles 
and resonances was discovered 0. Gell-Mann [Oj and Ne'eman [2j noticed that 
particles sharing the same quantum numbers (spin, parity) follow the symmetry of 
the mathematical group SU(3) which is based on 3 elementary generators, up, down, 
strange, or u, d, s, with spin 1/2 and fractional electrical charge, [1113 which Gell-Mann 
called quarks. Mesons are described as states made of a quark-anti-quark {qq) pair and 
baryons as states of 3 quarks (qqq) . This led to the prediction of a new baryon, the fl~ 
(sss) with strangeness -3, which was observed shortly thereafter. ITU] However, the Q~ 
had a problem: 3 identical s quarks in the same state, apparently violating the Pauli 
Exclusion Principle. To avoid this problem, it was proposed ^T] that quarks come in 
3 'colors', i.e. distinguishing characteristics which would allow 3 otherwise identical 
quarks to occupy the same state (formally, para- Fermi statistics of rank 3). A major 
breakthrough was the realization that the real SU(3) symmetry was not the original 
3 quarks uds (now called 'flavor'), but the 3 colors; and that color-charged gluons are 
the quanta of the 'asymtotically-free' strong interaction which binds hadrons. j2] The 
fourth "charm" or c quark, proposed to explain the absence of certain channels in 
weak decays of strange-particles, [12] thus had no problem fitting into this scheme-the 
quark symmetry became SU(2), groups of doublets. Elegant as these theories were, the 
experimental results were what made them believable. 

There were 3 key experimental observations that made the composite theory of 
hadrons believable: 1) the discovery of pointlike constituents ('partons') inside the 
proton, in deeply inelastic (large energy loss, z/, large 4-momentum transfer, Q) electron- 
proton (ep) scattering (DIS) at the Stanford Linear Accelerator Center (SLAC) [T3]; 2) 
the observation of particle production at large transverse momenta [pt) in p-p collisions 
at the CERN-Intersecting Storage Rings (ISR) [HI [121 [IHl, which proved that the 
partons of DIS interacted much more strongly with each other than the electromagnetic 
scattering observed at SLAC; 3) the observation of the J/\l/, a narrow bound state of cc, 
in both p+Be collisions at the BNL-AGS 12^, and in e~^e~ annihilations at SLAC |18j . 
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shortly followed by observation of the \E'', a similar state with higher mass ^H] which 
corresponded to cc bound states in a simple couloumb-like potential with a string-like 
linear confining potential. pUl 1^ |22] These discoveries turned Gell-Mann and Zweig's 
quarks from mere mathematical concepts to the fundamental constituents of matter, 
the components of the nucleon. 123 1211 1211 

1.1. From Bjorken Scaling to QCD to the QGP 

The fundamental idea to emerge from DIS, which was the basis of much of the subsequent 
theoretical developments leading to QCD, was the concept of Bjorken scaling |27, which 
indicated that protons consist of point-like objects (partons). The structure function 
-p2(Q^, ^) which describes the inelastic ep scattering cross section was predicted to |2Z1 
and observed to ^3] "scale" i.e. to be a function only of the ratio of the variables, 
Q'^/u. The deeply inelastic scattering of an electron from a proton is simply incoherent 
quasi-elastic scattering of the electron from point-like partons of effective mass Mx, 
where u = Q'^/2Mx, where M is the rest mass of the proton. Thus Bjorken 'x' is the 
fraction of the nucleon momentum (or mass) carried by the parton. Similar ideas for the 
scaling of longitudinal momentum distributions in p-p collisions were also given [^I29j. 
However these ideas related to the "soft" (low px) particle production rather than the 
large pt or "hard scattering" processes described by Bjorken [HO] . 

Bjorken scaling was the basis of QCD [2], the MIT Bag model jSI] and led to 
the conclusion [32] that "superdense matter (found in neutron-star cores, exploding 
black holes, and the early big-bang universe) consists of quarks rather than of hadrons" , 
because the hadrons overlap and their individuality is confused. Collins and Perry 
called this state "quark soup" but used the equation of state of a gas of free massless 
quarks from which the interacting gluons acquire an effective mass, which provides 
long-range screening. They anticipated superfluidity and superconductivity in nuclear 
matter at high densities and low temperatures. They also pointed out that for the 
theory of strong interactions (QCD), "high density matter is the second situation where 
one expects to be able to make reliable calculations — the first is Bjorken scaling" . In the 
Bjorken scaling region, the theory is asymptotically free at large momentum transfers 
while in high-density nuclear matter long range interactions are screened by many-body 
effects, so they can be ignored and short distance behavior can be calculated with the 
asymptotically-free QCD and relativistic many-body theory. Shuryak [HH] codified and 
elaborated on these ideas and provided the name "QCD (or quark-gluon) plasma" for 
"this phase of matter" , a plasma being an ionized gas. 

1.2. Relativistic Heavy Ions Collisions 

It was soon realized that the collisions of relativistic heavy ions could provide the means 
of obtaining superdense nuclear matter in the laboratory. [SH ESI EEl I^ZI The kinetic 
energy of the incident projectiles would be dissipated in the large volume of nuclear 
matter involved in the reaction. The system is expected to come to equilibrium, thus 
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heating and compressing the nuclear matter so that it undergoes a phase transition from 
a state of nucleons containing bound quarks and gluons to a state of deconfined quarks 
and gluons, the Quark Gluon Plasma, in chemical and thermal equilibrium, covering 
the entire volume of the colliding nuclei or a volume that corresponds to many units of 
the characteristic length scale. In the terminology of high energy physics, this is called 
a "soft" process, related to the QCD confinement scale [SHj 

AQ^D^(0.2GeV)-^^lfm . (1) 

Two energy regimes are discussed for the QGP At lower energies, typical of the 
AGS fixed target program, the colliding nuclei are expected to stop each other, leading 
to a baryon-rich system. This will be the region of maximum baryon density. At very 
high energy, 100 to 200 GeV per nucleon pair in the center of mass, nuclei become 
transparent and the nuclear fragments will be well separated from a central region of 
particle production. This is the region of the baryon-free or gluon plasma. In the nuclear 
fragmentation regions a baryon-rich plasma may also be formed jHIH EH] . 

There has been considerable work over the past three decades in making 
quantitative predictions for the QGP [S21 HH 1121 UHl The predicted transition 

temperature from a state of hadrons to the QGP varies, from ~ 150 MeV at zero 
baryon density, to zero temperature at a critical baryon density ^\ E5] roughly 1 
GeV/fm'^, ~ 6.5 times the normal density of cold nuclear matter, po = 0.14 nucleons /fm^, 
fiB — 930 MeV, where /zb is the Baryon chemical potential. A typical expected phase 
diagram of nuclear matter jlHI is shown in Fig. Not distinguished on Fig. ^ in the 
hadronic phase are the liquid self-bound ground state of nuclear matter and the gas of 
free nucleons ji7] . 
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Figure 1. a) (left) Proposed phase diagram for nuclear matter 0^1 : Temperature, T, 
vs Baryon Chemical Potential, /i. b) (right) Lattice calculation of energy density, 
e/T^ as a function of the number of active flavors: 2 flavor (u, d), 3 flavor (it, d, s). 

Predictions for the transition temperature for /i^ ~ are constrained to a relatively 
narrow range 140 < Tc < 250 MeV, while the critical energy density is predicted to be 
5 to 20 times the normal nuclear energy density, eo = 0.14 GeV/fm'^. Presumably, the 
most accurate predictions of the phase transition are given by numerical solutions of the 
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QCD Lagrangian on a lattice |1H1, see Fig. [TJd jl^l- Here, the critical energy density at 
the transition temperature, ~ 150 — 170 MeV, is stated to be presently known only 
with large errors jSU], ec=(0.3-1.3) GeV/fm^. 

One of the nice features of the search for the QGP is that it requires the integrated 
use of many disciplines in Physics: High Energy Particle Physics, Nuclear Physics, 
Relativistic Mechanics, Quantum Statistical Mechanics, etc. jHIl E21 EH] From the point 
of view of an experimentalist there are two major questions in this field. The first is 
how to relate the thermodynamical properties (temperature, energy density, entropy ...) 
of the QGP or hot nuclear matter to properties that can be measured in the laboratory. 
The second question is how the QGP can be detected. 

One of the major challenges in this field is to find signatures that are unique to the 
QGP so that this new state of matter can be distinguished from the "ordinary physics" 
of relativistic nuclear collisions. Another more general challenge is to find effects which 
are specific to A+A collisions, such as collective or coherent phenomena, in distinction to 
cases for which A+A collisions can be considered as merely an incoherent superposition 
of nucleon-nucleon collisions [Ml ISSl EE] • 

Many signatures of the Quark Gluon Plasma jSZl EH EH] have been proposed 
over the past two decades, which cover the experimental RHI programs at the AGS, 
the SPS, RHIC and soon the Large Hadron Collider (LHC) at CERN. In recent peer 
reviewed articles summarizing the first 3 years of RHIC operation [jHO I IF)T | IH ^ IH! ^ I(i4j . the 
four experiments and many distinguished theorists presented their results and opinions 
on QGP signatures and whether the QGP had been detected at RHIC. These were 
summarized in a BNL press release on April 18, 2005: "instead of behaving like a gas 
of free quarks or and gluons, as was expected, the matter created in RHIC's heavy 
ion collisions appears to be more like a liquid." This matter interacted much more 
strongly than expected, causing the theorists [HI] to give it the new name "sQGP" 
(strongly interacting QGP). Weighing heavily on this process was the CERN press 
release jHS] and unpublished preprint jHE] on February 10, 2000, just at the start of RHIC 
operations, which discussed properties of the QGP and announced that "The collected 
data from the [SPS] experiments gives compelling evidence that a new state of matter 
has been created. This state of matter found in heavy ion collisions at the SPS features 
many of the characteristics of the theoretically predicted quark-gluon plasma..." "The 
data from any one experiment is not enough to give the full picture but the combined 
results from all experiments agree and fit. Whereas all attempts to explain them using 
established particle interactions have failed, many of the observations are consistent 
with the predicted signatures of a quark-gluon plasma." In other words, several features 
of the CERN measurements were consistent with the expected properties of the QGP 
at that time, i.e. a gas of quarks and gluons. (Although not mentioned in any press 
release, the same could be said for several features of measurements at the AGS fixed 
target program.) In light of this apparent contradiction, the data from the AGS, SPS 
and RHIC experiments will be examined with emphasis on the QGP signatures outlined 
in these review articles and press releases. 
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Figure 2. a) (left) Schematic of collision of two nuclei with radius R and impact 
parameter b. b) (center) A p-p collision in the STAR detector viewed along the collision 
axis; c) (right) Au+Au central collision at ^/snn = 200 GeV in the STAR detector. 



A schematic drawing of a Relativistic Heavy Ion collision is shown in Fig. |2t^. In 
the center of mass system of the nucleus-nucleus collision, the two Lorentz-contracted 
nuclei of radius R approach each other with impact parameter b. In the region of 
overlap, the "participating" nucleons interact with each other, while in the non-overlap 
region, the "spectator" nucleons simply continue on their original trajectories and can 
be measured in Zero Degree Calorimeters (ZDC), so that the number of participants can 
be determined. The degree of overlap is called the centrality of the collision, with 6 ~ 0, 
being the most central and b ~ 2R, the most peripheral. The maximum time of overlap 
is To = 2/2/7 c where 7 is the Lorentz factor and c is the velocity of light. The energy 
of the inelastic collision is predominantly dissipated by multiple particle production. 

For any observed particle of momentum p, energy E, the momentum can be resolved 
into transverse (pt) and longitudinal (pi) components; and in many cases the mass (m) 
of the particle can be determined. The longitudinal momentum is conveniently expressed 
in terms of the rapidity (y): 

y - in {^) (2) 



coshy = E /rriT sinhy = Phlmx dy = dpi/ E (3) 

where 

rriT = \J + p^ and E = \Jp\ + rn^ = \Jp'^ + rn^ (4) 
In the limit when (m ^ E) the rapidity reduces to the pseudorapidity {rf) 

77 = -lntan^/2 (5) 

cosh 7] = CSC 6 sinh 77 = cot 6 (6) 



where 6 is the polar angle of emission. The rapidity variable has the useful property 
that it is additive under a Lorentz transformation. 
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2.1. Kinematics of the collision 

For any collision, the center-of-mass (cm.) system — in which the momenta of the 
incident projectile and target are equal and opposite — is at rapidity y^^. The total 
energy in the cm. system is denoted a/s, which, evidently, is also the "invariant mass" 
of the cm. system. For a collision of an incident projectile of energy mass mi, 
in the "Laboratory System", where the target, of mass m2, is at rest (appropriate for 
fixed-target experiments): 

s = ml + + 2Eim2 . (7) 

The cm. rest frame moves in the laboratory system (along the collision axis) with a 
velocity jS'^^c corresponding to: 

^cm _ El +1712 ^cm _ ^^^^-1 ^cm ^ ^g) 

where 7^ = 1/(1 — /3^). Another useful quantity is y^'^^'^^ the rapidity of the incident 
particle in the laboratory system 

nil 

and note that for equal mass projectile and target: 

^cm _ ^beam/2 . (10) 

In the region near the projectile or target rapidity, the Feynman x fragmentation variable 
is also used: 

XF = 2pl/V^ , (11) 

where is the longitudinal momentum of a particle in the cm frame. 

The kinematics are considerably simpler in the nucleon-nucleon cm. frame in 
which the two nuclei approach each other with the same 7. (This is the reference 
frame of the detectors at RHIC). The nucleon-nucleon cm. energy is denoted ^/snn, 
and the total cm. energy is ^/s = A ■ ^/snn for symmetric A+A collisions. The 
colliding nucleons approach each other with energy ■s/snn/2 and equal and opposite 
momenta. The rapidity of the nucleon-nucleon center of mass is = 0, and, taking 
mi = m2 = mAT =931 MeV, the projectile and target nucleons are at equal and opposite 
rapidities 



^proj _ _^target _ ^^^j^-l VfN^ ^ . (12) 

2m]y 

2.2. A brief overview of relativistic collisions of nucleons and nuclei 

The challenge of RHI collisions can be understood from Fig. Eb, which compares a p-p 
collision to an Au+Au central collision in the STAR detector jH^]. It would appear to 
be a daunting task to reconstruct all the particles in such events. Consequently, it is 
more common to use single-particle or multi-particle inclusive variables to analyze these 
reactions. 
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A single particle "inclusive" reaction involves the measurement of just one particle 
coming out of a reaction, 

a + 6 — > c + anything . 

The terminology j2Hl comes from the fact that all final states with the particle c are 
summed over, or included. A "semi-inclusive" reactionjnZI refers to the measurement of 
all events of a given topology or class, e.g. 

a + b ^ 111 particles in class 1 + anything , 

where "centrality" is the most common class in relativistic heavy ion collisions. 

Measurements are presented in terms of the (Lorentz) invariant single particle 
inclusive differential cross section (or Yield per event in the class if semi-inclusive): 

Ed^a d^a 1 



(13) 



dp^ PTdprdydcf) 2n 
where y is the rapidity, pt is the transverse momentum, and (j) is the azimuth of 
the particle (see Fig. The average transverse momentum, (pr), or the mean 
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Figure 3. a) (left) Semi-inclusive invariant pT spectra for tt*, K^, in Au+Au 
collisions at ^/snn — 200 GeV [HSI- b) (right) (pr) of positive particles as a function 
of centrality (A^part) from the same data. 



transverse kinetic energy, (m-r) — m, or the asymptotic slope are taken as measures 
of the temperature, T, of the reaction. 

It is important to be aware that the integral of the single particle inclusive cross 
section over all the variables is not equal to aj the interaction cross section, but rather 
is equal to the mean multiplicity times the interaction cross section: (n) x aj. Hence 
the mean multiplicity per interaction is 



(n) 



1 



2n 



dydpT Pt i{PT,y) 



1 



dy 



da 
dy 



dyp{y) 



(14) 
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where the terminology for the muhiphcity density in rapidity is (I/ct/) da/dy = p{y) = 
dn/dy for identified particles [m known), dn/drj for non-identified particles (m unknown, 
assumed massless). The total charged particle multiplicity is taken as a measure of the 
total entropy, S and dn/drj is taken as a measure of the entropy density in restricted 
intervals of rapidity. 



2.2.1. The rapidity density and the rapidity plateau. The shape and evolution with a/s 
of the charged particle density in rapidity, dn/dy, provide a graphic description of high 
energy collisions. Data from a classical measurement in a streamer chamber from p-p 
collisions at the CERN ISR are shown in Fig. HJ Regions of nuclear fragmentation 
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Figure 4. Measurements in a streamer chamber at the CERN ISR of the 
normalized charged particle densities -^^^ (correctedfor acceptance up to \rj\ ~ 4) in 
various intervals of the total observed multiplicity as a function of the cm. energy ^/s 
of the p-p collision. 



take up the first 1-2 units around the projectile and target rapidity and if the center- 
of-mass energy is sufficiently high, a central plateau is exhibited. The data in p+A and 
A+A collisions follow a similar trend (Fig|3^) [70j. The distributions increase in width 
with increasing ^/snn but by a smaller amount than the increase in y^^^™ and show 
a small decrease in width with increasing centrality. dn/drj increases with increasing 
centrality, ^/snn and A in A+A collisions. In the asymmetric d+Au collision, dn/drj in 
the target rapidity region is larger than in the projectile region, but not by much, only 
about 50%. Also the nuclear transparency is evident, there is no reduction of particles 
at the projectile rapidity with increasing centrality. 
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Figure 5. a) (top right) dn/drj for A+A and d+Au collisions at RHIC as a 
function of ^/snn IZO]- b) (bottom left) dn/dy\p — dn/dy\p at AGS, SPS and RHIC, 
yP™j = 1.6,2.9,5.4 [HI- 



Subtleties of the distributions in A+A collisions become apparent when identified 
particles are used [7T]. In Fig. the difference of dn/dy for protons and anti-protons, 
i.e net-protons is shown as a function of cm. energy, yijvTv = 5 (AGS, Au+Au), 17 
(SPS, Pb+Pb), 200 (RHIC, Au+Au), y^'"'^ =1.6, 2.9, 5.4. As is reduced, stopping 

of the participating nucleons is indicated by the nuclear fragmentation peak moving from 
the fragmentation region (not visible for RHIC) to mid-rapidity. 

2.2.2. The Bjorken energy density Another variable, closely related to multiplicity, is 
the transverse energy density in rapidity or dEx/dy ~ {px) x dn/dy, usually measured 
in calorimeters by summing over all particles on an event in a fixed but relatively large 
solid angle [12]: Et = J2iEi sin 9i. dEx/dy is thought to be related to the co- moving 
energy density in a longitudinal expansion [731 113 , and taken by all experiments as a 
measure of the energy density in space e: 

- ^^^^ 
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where Tp, the formation time, is usually taken as 1 fm/c, ttR'^ is the effective area of 
the collision, and d (Et) /dy is the co-moving energy density. 



E866 & E877-AGS corr to At} = 1 and A*=27T 14.6 A GeV/c PHENIX and EB02 E, Transverse Energy corr to A7? = l and A*=27r 




E, (A*.27,) (GeV) (GeV) 

Figure 6. a) (left) Mid-rapidity Et spectra (corrected to A0 = 27r, Arj — 1) from AGS 
measurements by E814/E877 ■ and E802/E866 • O j^g] for 14.6-^ GeV/c Si+Al, 
Au and 11.6-A GeV/c Au+Au (corrected to 14.6-^ GeV/c). The E802 data have been 
scaled by the factor indicated to match the E877 measurement. Upper percentiles of 
the distribution are indicated. The light O points on the Au+Au data are from a ZDC 
trigger of (0-7%) centrality. b) (right) Et distribution for y^sj^ = 200 GeV Au+Au 
at RHIC • [77] together with the E802 measurement O from the left panel scaled 
as indicated |Z8|. The solid circles represent a measurement in an azimuthal aperture 
of A(j) = 5 X 22.4°, corrected to full azimuth. The other solid symbols represent 
measurements in A(p = 4, 3,2, 1 x 22.4°, corrected to full azimuth, where the smaller 
the aperture, the flatter the slope of the data above the knee. 

Besides the Bjorken energy density, the importance of Et distributions in RHI 
collisions (see Fig. ^ is that they are sensitive primarily to the nuclear geometry of the 
reaction, and hence can be used to measure the centrality of individual interactions on 
an event-by-event basis [71j. Fig. IHt- shows mid-rapidity Et distributions for Si-|-Al, 
Si-|-Au and Au-|-Au at the AGS [721 IZB] • The increase of dET / drj with increasing atomic 
mass of projectile and target is evident. Fig. [Hb shows that Et distributions at RHIC 
and the AGS are the same shape, when scaled to match at the knee, which shows that 
the shape of the distribution is essentially independent of ^/snn and thus dominated by 
the nuclear geometry of the Au-|-Au reaction. Above the knee, the slope of the fall-off 
is sensitive to fluctuations and depends somewhat on the solid angle [7^1 IZ3 ESI- ^ 
consistent evaluation [77| [7H1 of the Bjorken energy density for 0-5% centrality gives 
TpeBj = 1.0, 2.0, 5.4 GeV/fm2 at = 5, 17, 200 GeV. 

2.3. Space-time and quantum mechanical issues 

An interaction of two relativistic heavy ions can be viewed initially as the superposition 
of successive collisions of the participating nucleons in the overlap region (recall 
Fig. 121) ffS^. Collective effects may subsequently develop due to rescattering of the 
participants with each other or with produced particles. Conceivably, a cascade of 
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interacting particles could develop. However, the actual situation is considerably more 
fascinating. 

Space-time and quantum mechanical issues play an important role in the physics 
of RHI collisions and the considerations are different in the longitudinal and transverse 
directions. In QED and QCD, which have the same 1/r form of the potential at small 
distances, two particles of charge Zi and Z2, velocity v, which pass each other at impact 
parameter b each acquire a momentum pt = 2ZiZ2/bv transverse to the direction of 
motion. Thus large transverse momenta correspond to small impact parameters in both 
QED and QCD. 

In the longitudinal direction because of the large 7 factors involved, small 
excitations/deexcitations of the colliding nucleons, take place over long distances. When 
a nucleon with momentum, pi, mass, M, makes a collision, the only thing it can do 
consistent with relativity and quantum mechanics is to get excited to a state with 
invariant mass M* > M, with roughly the same energy and reduced p'l^ = Pl — ^Pl, 
where Ap| = — Am|. ~ — Am^ from Eq. 01 By the uncertainty, principle, a distance 
6z = h/Api = 7/5/Am(=14 fm for 7 = 10 and m = m.,^) is required to resolve Am. The 
large 7 factor in relativistic collisions ensures that the excited nucleons pass through the 
entire target nucleus before de-exciting into e.g. a nucleon + a pion. Thus, nuclei are 
transparent to relativistic nucleons; and pions are produced outside the target nucleus, 
thus avoiding a nuclear cascade. 

For instance, in the collision of a relativistic proton with a 15 interaction-length- 
thick lead brick, a cascade develops and all particles are absorbed. Nothing comes out 
the back. By contrast, in the collision of a relativistic proton with a lead nucleus, which 
is roughly 15 interaction mean-free-paths thick through the center, the (excited) proton 
comes out the back! This is relativity and quantum mechanics in action. 



2.4- Participant fTVpartj scaling-the Wounded Nucleon Model 

These concepts of relativity and quantum mechanics are dramatically illustrated by 
particle production in proton nucleus (p+A) interactions measured at bombarding 
energy of 200 GeV ^/SNN = 19.4 GeV. When a high energy proton (or any hadron) 
passes through a nucleus, it can make several successive collisions. However, the 
charged particle multiplicity density, dn/dy, observed in p+A interactions is not simply 
proportional to the number of collisions, but increases much more slowly. These features 
are strikingly illustrated jHOj in Fig- HI which is the pseudorapidity distribution of 
relativistic charged particles {v/c > 0.85) produced by 200 GeV proton interactions in 
various nuclear targets. The sizes of the nuclei are discussed in terms of z/, the average 
number of collisions encountered by an incident hadron passing through a nucleus of 
atomic mass A in which an interaction occurred: 

z/ = , (16) 

where ahp and ahA are the absorption cross sections for the incident hadron on a nucleon 
and a nucleus, respectively 
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Figure 7. a) (left) Pseudorapidity distributions dn/drj of relativistic charged particles 
for various values of v in 200 GeV/c proton- nucleus interactions [HD]- b) (center) Ra 
as a function of v for collisions of 200 GeV tt^ on various nuclei. The solid line 
shows the linear trend of the data ^ Eq. El !^ • c) (right) Mean number of pions per 
participant for A+A central collisions and p{p) — p collisions as a function of ^/snn 
represented by F = (^s^ - 2mjv)3/4/^i]^i/4 [gg 



The most dramatic feature of Fig. [7| is that there is virtually no change in the 
forward fragmentation region [r] > 5.0) with increasing A. By contrast, there is 
tremendous activity in the target region (r/ < 0.5). In the mid-rapidity region, dn/drj 
increases with A and the peak of the distribution shifts backwards. The integral of the 
distribution (the average multiplicity) and dn/drj at mid-rapidity show a linear increase 
with v for all targets given by the relation 

which is simply the ratio of the average number of participating nucleons (participants 
or A/'part) for the 2 cases: A'part = 1 -|- i> for the p-|-A collision and A^'part = 2 for a 
p-|-p collision. In addition to the quantum and relativistic arguments given above, this 
result can be explained by the further assumptions that the excited nucleon interacts 
with the same cross section as an unexcited nucleon and that the successive collisions of 
the excited nucleon do not greatly affect the excited state or its eventual fragmentation 
products. jHSl This leads to the conclusion that the multiplicity in nuclear interactions 
should be proportional to the total number of projectile and target participants, rather 
than to the total number of collisions. This is called the Wounded Nucleon Model 
(WNM) [Bl. 

Interestingly, the WNM works well only at roughly ^/sNlv ~ 20 GeV where it was 
discovered. At lower ^/snn, at the AGS, particle production in p-|-A and A-l-A collisions 
is smaller than the WNM and for > 31 GeV is larger than the WNM jHSl IHEl EH 
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as emphasized by Gazdzicki [H21 ( "the kink"). However as the effect is seen in p+A 
colhsions, it is not hkely to be related to the QGP [SHI 1111 HDI although the physics is 
surely very interesting and needs further investigation. 

3. Signatures of the Quark-Gluon Plasma 

One of the more interesting signatures proposed for the QGP is that it could trigger a 
catastrophic transition from the metastable vacuum of the present universe to a lower 
energy state [23121], "a possibility naturally occurring in many spontaneously broken 
quantum field theories" . This type signature gives laboratory directors nightmares, and 
such a possibility must be ruled out to a high degree of certainty |S2]- However this 
is actually one example of a class of signatures of the QGP called "Chiral Symmetry 
Restoration" . It is convenient to group QGP signatures into classes according to their 
sensitivity to expected properties of the QGP [SHI: (i) signals sensitive to the equation 
of state; (ii) signals of chiral symmetry restoration; (iii) probes of the response function 
of the medium (including deconfinement). 

All authors imagine that the QGP is in a state of chemical and thermal equilibrium 
at temperature Tp which is established over a short formation time, Tp, following the 
collision and that hydrodynamics describes the evolution of the system of quarks and 
gluons as it expands and cools and eventually freezes out at transition temperature, 
Tc, at time, Tc, to a system of hadrons which are observable. Another important time 
is the overlap time of the colliding Lorentz contracted nuclear pancakes as seen in the 
nucleus-nucleus cm. system, tq = 2R/j, where R is the radius of the nucleus in a 
symmetric A+A collision. For Au+Au or Pb+Pb collisions, R ~ 7.5 fm , ro = 5.6, 1.7, 
0.14 fm/c at =5, 17, 200 GeV. 

3.1. Degrees of freedom, equation of state 

It is important to emphasize that the operative degrees of freedom in the QGP are 
taken to be quarks and gluons and the operative "charge" is taken be color, unscreened 
over many units of the characteristic length scale. Since the spin-1 gluons have 8 colors 
and the spin-1/2 quarks have 3 colors, with presumably 3 active flavors in the QGP 
{u,d,s), there are many more degrees of freedom in the QGP than in the final state 
'gas' of hadrons (HG); and this is one of the signatures of the transition from one 
phase to the other. The QGP is described by the thermodynamic quantities, entropy 
density, s, energy density, e, temperature, T, pressure P, volume, V, and it is important 
to understand how to relate these quantities to properties of the observable final state 
particles. As in classical thermodynamics, the equation of state (EOS) is the relationship 
among the three quantities, P,V,T, which are not independent for a fixed amount of 
matter. A relation between any two quantities is sufficient to define the EOS. For 
instance, the large difference in degrees of freedom [63^ between a QGP (taken as a 
Stefan-Boltzmann gas) and a hadron gas dominated by spin pions (3 d.o.f: tt^, tt^, tt") 
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is shown in Lattice calculations (Fig. [TJd) 



30 30 

However, it is not obvious how this effect could be measured — possibly by evidence of 
discontinuities or rapid changes as a function of an observable jOHl IHll or by fluctuations 
due to passing through the phase transition j^Eni- (Recall that T ^ {pt), s dn/drj 
and e ^ dET/dy). 

In general, fluctuations are a standard probe of thermodynamic systems and phase 
transitions [HZl IHHl UHl HHI and have been extensively studied at the SPS and RHIC. 
However, no strong signals have yet been observed in relativistic A+A collisions. |l()()j . 

The volume of a thermalized source is thought to be measurable by identical 
particle interferometry using the GGLP effect |101j . commonly called Bose- Einstein 
correlations since the measurement is predominantly performed with identical bosons 
{7r^,K^) |102[ llOtjj . Also, collective effects, where all partons (or particles) moving in 
the same direction have a common (flow) velocity, are sensitive to the EOS. |1()4[ \U)^ 

3.1.1. Thermal Equilibrium One of the best probes of thermal equilibrium is thermal 
lepton pair production [lOH] or thermal photon production |107j . This is effectively the 
'black-body radiation' of the QGP and should follow a Boltzmann distribution in the 
local rest frame 



"^'^ ^ ~e-^/^ . (19) 



dpiPrdpT e-^/'^ ± 1 

Since prdpT = rriTdmT and E = rriT coshy, a signal of thermal production is that 
the Pt and mass dependence of the cross section are not independent but depend 
only on the transverse mass m-p. This means that at any fixed value of and 
rapidity, the {pr^iTiee)) is linearly proportional to ruee- As neither the photons nor e+e" 
pairs are strongly interacting, they emerge from a QGP or hadronic system without 
interacting and thus are sensitive to the entire thermal history of the system, especially 
the early stage where the QGP should be dominant. In principle, the initial temperature 
of the system, T, can be determined from the rate of thermal photon or e^e~ pair 
production jl()7^ llOf)] . However, even though radiation from the QGP comes from the 
interaction of 'massless' quarks and gluons: 

gg— >e^e^ g q ^ 1 q qq^qq'j etc, 

the latest 'state-of-the-art' calculations |109[ IllOj indicate that for the same 
temperature, the thermal photon production rates in a QGP or a hadron gas are 
very similar |lllj . It is important to distinguish thermal photons (lepton pairs), which 
dominate at low px (mass), near Tc, and are exponential, from the high px (mass) direct 
photons (Drell-Yan pairs) produced by hard-scattering in some of the same subprocesses, 
Qg^lQiQQ^ c^e~), which follow a power-law. Hard-scattering (see below) is not 
relevant for thermal equilibrium. 
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3.2. Chiral Symmetry Restoration 

In the QGP |112[ I3H] the hght quarks are expected to have the same mass or be massless, 
thus exhibiting exact "chiral symmetry". These are the so-called "current quarks" of 
DIS. This is quite different from the normal vacuum state of QCD in which chiral 
symmetry is "spontaneously broken" and the quarks are massive, called "constituent 
quarks" . This leads to several interesting effects which might be measurable. 

3.2.1. Strangeness Enhancement In the QGP, the gluons, quarks and anti-quarks 
continuously react with each other via the QCD subprocesses: 

gg^qq qq^gg qq^q'q' 

where q' represents a different flavor quark {u,doTs). After several interactions have 
taken place, the reaction rates and the abundances of the gluons and the different flavor 
quarks (and anti-quarks) will become equilibrated, so that they no longer change with 
time. This is called chemical equilibrium. Since the masses of u, d, s quarks should 
be nearly degenerate if chiral symmetry is restored, they should reach nearly the same 
equilibrium value so that the strange quarks s, s should have nearly the same abundance 
as the u, it and d, d in the gluon plasma. Thus strangeness should be enhanced in the 
QGP |llH^lll4(lll5j compared to p-p collisions where the abundance of strange particles, 
e.g. K^, is much below their thermal equilibrium value relative to vr^. In the baryon- 
rich plasma, the s, s will be enhanced compared to u and d, since u, u and d, d are 
"Pauli" blocked by valence u and d quarks. 

Since the quarks and gluons in a QGP can not be observed directly, the principal 
probe of chemical equilibrium and strangeness enhancement in the QGP is the particle 
composition of observed hadrons. For instance, the abundance of strange mesons and 
baryons as well as anti-baryons should be quite different in a QGP than in a hadron gas 
or in an ordinary nuclear collision. 

3.2.2. Disordered Chiral Condensate An interesting anomaly possible when the 
medium returns to the normal QCD ground state from the chirally symmetric QGP state 
is that "misaligned" QCD vacuum regions might occur such that instead of emitting 
7r° and it~^tt~ with a binomial distribution with probability of 1/3 for the fraction 
R = 7r°/(7r° + 7r+ + tt^) of 7r° emitted, individual "misahgned" domains could emit, 
for instance, only tt^ or only vt+tt^ pairs |116j . This would give a distribution of the 
form P{R) = 1/(2/R), which has the same mean value, (R) = 1/3 but a totally different 
event-by-event distribution. Incredibly, so far there is only one (null) measurement |117j 
(at the SPS) due to the inherent difficulty of reconstructing 7r° — > 7 + 7 at low pt- 

3.2.3. mass-shifts, branching ratio changes While it is clear how chiral symmetry 
breaking and restoration affect the quarks, it is not so clear how they affect the particles 
we observe which are by definition on their mass-shell, with their well defined masses and 
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are probably formed during the hadronization process, well after the QGP has cooled. 
Pisarski |118j has suggested that if the temperature for chiral symmetry restoration, T^h 
is greater that the QGP critical temperature, T^, then all hadrons except pions would 
unbind at T > Tc, with pions massive quarks and gluons in coexistence until T > Tch, 
when all quarks and gluons become massless. This might mean a state of constituent 
quarks in the region Tc < T < Tch- Also, particles which are formed and which decay to 
leptons or photons in this phase can be detected, since the leptons and photons do not 
interact with the QGP and simply emerge unscathed. As the width of the p is 150 MeV 
(lifetime 1.3 fm/c), this might give rise to |118j an "extraordinary signal" of a thermal 
p° — s> dilepton peak "quite distinct from the familiar p° peak". To a lesser extent, a 
similar effect is predicted for the (p meson |119j which has a width of 4.3 MeV (lifetime 
~ 46 fm/c). The mass is just above the K^K~ mass but its width is dominated by 
K^K~ (or KiKg) decay and therefore is very sensitive to relative shifts of its mass 
with respect to the masses of the Kaons. The width of the is predicted to increase by 
a factor of 2-3 due to this effect, which may be detectable directly by high resolution 
measurements of — > dileptons, or by an increase of the observed (p ~^ e^e~ /K^K~ 
branching ratio from the standard (vacuum) value. 

3.3. Probes of the response function of the medium including deconfinement 

3.3.1. Deconfinement Since 1986, the 'gold-plated' signature of deconfinement was 
thought to be J/\l/ suppression. Matsui and Satz proposed |12Uj that J/\E' production 
in A+A collisions will be suppressed by Debye screening of the quark color charge in the 
QGP. The J/\E' is produced when two gluons interact to produce a c, c pair which then 
resonates to form the J /'^ . In the plasma the c, c interaction is screened so that the c, c 
go their separate ways and eventually pick up other quarks at the periphery to become 
open charm. Due to the fact that it takes time for the initial c, c pair to form a J/^', 
for the quarks to separate to the correct Bohr orbit, the J /'^ suppression should vanish 
characteristically with increasing px as the c, c pair leaves the medium before the J/^f 
is formed or screened jl21j . J/^^ suppression would be quite a spectacular effect since 
the naive expectation was that J /'^ production, due to the relatively large ~ 1.5 GeV 
scale of the charm quark mass, should behave like a pointlike process, proportional to 
y4 X A in an A+A collision, and thus would be enhanced relative to the total interaction 
cross section, which increases only as A^^'^ . 

The screening of the coulomb-like QCD potential for heavy quarks is supported by 
lattice gauge calculations |122| l3n] (Fig. |H1). With increasing temperature, T, in analogy 
to increasing Q^, the strong coupling constant agiT) becomes smaller, reducing the 
binding energy, and the string tension, cr(T), becomes smaller, increasing the confining 
radius, effectively screening the potential |I3]: 

V(r) = J-"-^ + ar^-'-"-±e->" + a^l^^ (20) 

3 r 3 r p 

where p = p(T) = l/ro'is the Debye screening mass. jlS]- When the potential in Fig.|S^ 
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Figure 8. a) (left) Lattice gauge calculations of heavy quark potential V{r) as a 
function of T/T^ The solid line is the normal (T = 0) potential; b) (right) 

Estimate |50| of rmed, the distance beyond which the force between a static quark 
anti-quark pair is strongly modified by temperature effects, compared to the Debye 
screening radius, rn — and the mean square charge radii of cc and bb bound 

states (lines). Note that the actual melting points of the 4* , Xc and J/^ are no longer 
expected to be ^0.2,0.7 and 1.1 Tc as in the calculation illustrated [^01, but rather 1.1, 
1.1 and 2Tc by more recent calculations [HSl El ESI • 



becomes constant with increasing radius, the binding force vanishes as illustrated [SHj 
in Fig. IHb- 

Since the radii of higher excitations of the bound cc (charmonium) and bb states 
increase with their masses and the binding energies decrease, the states should melt 
sequentially with increasing T from the highest mass to the lowest in a given family 
(cc, bb). Since 60% of the observed J/\I' are directly produced, 30% are from Xc 
decay and 10% from ^' decay, this should give a characteristic sequential nature 
to J/^ suppression with increasing centrality (presumably increasing T) of an A+A 
collision — first the 10% of the J/'^ from the ^' are suppressed, then the 30% from 
the Xc are suppressed and eventually the other 60% vanish when the direct J/^ 
melts. Such an effect was apparently observed at the SPS [ISH] and wonderful models 
were woven to explain the results |127t I128j . However, recent measurements from 
RHIC (see below), and an increase in the expected dissociation temperatures for the 
charmonium states jl23U124[ I125j put the whole concept of J/\E' suppression as a probe 
of deconfinement into question. 

3.4- Jet Quenching 

A new tool for probing the color response function of the medium was developed in the 
early 1990's [T2ni IMIl ESI] and given a firm basis in QCD |in21 CHSl CBll CHSl CHSl just 
before RHIC turned on. In the initial collision of the two nuclei in an A+A collision 
when the Lorentz contracted nuclei are overlapped, hard scatterings can occur which 
produce pairs of high pt outgoing partons with their color charge fully exposed as they 
travel through the medium before they fragment to jets of particles. If the medium 
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also has a large density of color charges, then the partons are predicted to lose energy 
by 'coherent' (LPM) gluon bremsstrahlung which is sensitive to the properties of the 
medium. This leads to a reduction in the pt of both the partons and their fragments 
and hence a reduction in the number of partons or fragments at a given pt, which is 
called jet quenching. The effect is absent in p+A or d+A collisions due to the lack of a 
medium produced (see Fig. ^ 
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Figure 9. a) (left) Schematic diagram of liard scattering in A+A and d(p)+A 
collision. For Au+Au, high px partons from scattering when the nuclei overlap, 
emerge sideways through the medium formed. For d+Au no medium is formed and 
the outgoing partons travel in the vacuum until they fragment, b) (right) Transport 
coefficient g as a function of energy density for different media |136| : cold nuclear 
matter (■), massless hot pion gas (dotted) and 'ideal' QGP (solid curve) 



The screened coulomb potential f Eg. 1201) resulting from the thermal mass /i acquired 
by gluons in the medium modifies the elastic scattering of outgoing partons such that 
the average 4-momentum transfer per collision is finite: 



da, 



el 



dq" 



dq' 



(21) 



The integral in the denominator is aei = 1/pAmfp, where p is the density of scatterers 
and Amfp is the mean free path for elastic scattering. The 'transport coefficient' q is 
the mean 4-momentum transfer /collision expressed as mean 4-momentum transfer per 
mean free path, so that the mean 4-momentum transfer for length L in the medium 
is {q{LY) = qL = /x^L/Amfp. Any gluon radiation that is contained within the cone 
defined by the transverse momentum, (g(L)), is coherent over all the scatterings (LPM 
effect) so that the energy loss of a parton due to gluon bremsstrahlung per unit length 
(x) of the medium takes the form [136J: 

—^c:^as{q{Lf) = asqL = asP^L/\rni^ . (22) 
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Thus the total energy loss in the medium goes like |134j . It is also important to note 
that the screening mass, /i^, plays the role of tmin = Ql ^^e gluon bremsstrahlung in 
the medium |132j . 

The solution for the energy loss in pQCD for a static medium can be calculated 
analytically in only two limits, many soft scatterings (like multiple coulomb scattering) 
and 1 hard scattering, so that many treatments have been developed in varying degrees 
of approximation in the opacity, L/Amfp (hke u, Eg. IT^. the average number of partonic 
collisions in the medium |137| \1'AH\ lllfflj . Also, a formulation convenient for numerical 
simulations has been given |140j . However, the medium is not static but is expanding, 
flowing, etc, which gives rise to further complications. It is generally agreed that the 
energy loss measures q, which depends on the color charge density in the medium within 
the screening radius around the probe (Fig. Eb)- However it is not sensitive to the 
vanishing of the string tension beyond ro = which represents true 'extended' 

deconfinement because the formation time to resolve a pion in the fragmentation 
process means that quarks do not fragment to on-the-mass-shell pions until they are 
well outside medium where ordinary confinement applies. 

4. The Search for the QGP 

In the early 1980's, before measurements with relativistic heavy ions of cm. energy 
^/sNN ^ 2 GeV (~ 2 nucleon masses) became available, it was expected that the search 
for the QGP would proceed by looking for its predicted properties, or by the discovery 
of "anomalies" or discontinuities as a function of some experimental observable in A+A 
collisions jl41j . However, it was also realized that "systematic studies and comparison 
of pp, pA and AA data are equally important to understand basic processes hiding 
behind the phenomena observed in AA collisions" |142j . Thus it was not unreasonable 
to expect a few surprises: 

• Many of the predicted properties will be found, but will not be the QGP 

• The QGP will have a few unpredicted or unexpected properties 

• The search will uncover many unexpected backgrounds and new properties of pp 
and AA collisions, some of which may be very interesting pheonmena in their own 
right. 

Consequently, before proclaiming a perceived anomaly in A+A collisions as the 
QGP, it is imperative to learn and understand the "ordinary physics" of relativistic 
nuclear collisions as well as p-p and p+A collisions. As emphasized by Van Hove |141j 
much could be gained in this regard by studying the history of strong interactions for 
the past 30 years. An understanding of the properties of high energy p-p and p+A 
interactions is vital to the ability to distinguish the "ordinary physics" of relativistic 
nuclear interactions from the signatures of production of a new phenomenon like the 
QGP. It is also possible that the "ordinary physics" may in itself be quite interesting. 
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Furthermore, by studying the hard won knowedge of the past, one might hope to avoid 
some pitfalls in the future. 

4.1. p-p physics 

The theory and phenomenology of secondary particle production in ultrarelativistic 
hadron-hadron collisions originated with the study of cosmic rays |143j . It is hard to 
overstate the importance of the fundamental observation made by cosmic ray physicists, 
that the average transverse momentum of secondary particles is limited to ~0.5 GeV/c, 
independent of the primary energy |144t I143j . Cocconi, Koester and Perkins |145j 
proposed the prescient empirical formula for the transverse momentum spectrum of 
meson production: 



where px is the transverse momentum in GeV/c and (px) = 2/6 =0.333 GeV/c. The 
observation by Orear |146j that large angle p-p elastic scattering at AGS energies (10 to 
30 GeV incident energy) "can be fit by a single exponential in transverse momentum, 
and that this exponential is the very same exponential that describes the transverse 
momentum distribution of pions produced in nucleon-nucleon collisions", led to the 
interpretation |148j that particle production was "statistical", with Eq. |2Slas a thermal 
Boltzmann spectrum, with 1/6=0.167 GeV/c representing the "temperature", T at 
which the mesons or protons were emitted |149j . 

It was natural in a thermal scenario |150| I151j to represent the invariant cross 
section, as a function of the longitudinal variable rapidity, (y) , in terms of the transverse 



mass, m-r = ypj^ + m^, with a universal temperature parameter T. This description 
nicely explained the observed successively increasing (px) of vr, K, p, with increasing 
rest mass (e.g. see Fig. ^ and had the added advantage of explaining, by the simple 
factor Q-^i^K-rriTr) ^ 11% ^ the low value of ~10% observed for the K/ir ratio at low px 
at ISR energies (v^ ~ 20 - 60 GeV). 

The introduction of the constituent quark model, which used SU(3) symmetry 
to explain the hadron flavor spectrum and the static properties of hadrons jHl IH], as 
a dynamical model to calculate the flavor dependence of identified hadrons in soft 
multiparticle production, |152j together with the inclusive reaction formalism j2Hl 1211 
I153j . which showed that there was much to be learned by simply measuring a single 
particle spectrum, brought the study of identified inclusive single particle production 
into the mainstream of p-p physics. However, in the constituent quark model, the 
"suppression" of strange quarks, evident from the small K/tt ratio, was not explained, 
but simply quantified |154j by a parameter A = 2ss/{uu + dd), which represents the 
ratio of the numbers of ss pairs produced to the average of uu and dd pairs. 



da 



(23) 



PxdpT 




4- 1.1. Thermodynamics and Hydrodynamics in p-p collisions One of the burning issues 
in the early 1950 's in soft multiparticle production was whether more than one meson 
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could be produced in a single nucleon-nucleon collision ("multiple production") or 
whether the multiple meson production observed in nucleon-nucleus interactions was 
the result of several successive nucleon-nucleon collisions, with each collision producing 
only a single meson ("plural production") |155j . The issue was decided when multiple 
meson production was first observed in 1954 in collisions between neutrons with energies 
up to 2.2 GeV, produced at the Brookhaven Cosmotron, and protons in a hydrogen filled 
cloud chamber |156t I157j . 

The observation of multiparticle production occurring not only in nucleon-nucleus, 
but also in nucleon-nucleon collisions motivated Fermi and Landau to develop the 
statistical |158j and hydrodynamical |159j approach to multi-particle production. 
Belenkij and Landau observed that although the statistical model of Fermi is sufficient 
to describe the particle numbers in terms of only a temperature and a chemical potential, 
this model has to be extended to hydrodynamics, when particle spectra are also 
considered. To quote Carruthers |16Uj . "One envisions a thin slab of hot hadronic 
matter in thermal equilibrium just after the collision; strictly speaking this is a 'head- 
on' collision picture, but one can imagine that a fraction of the collision products are 
described by this initial condition while leading particles carry away a sizeable fraction 
(of the order of 1/2) of the energy and perhaps most of the angular momentum. In 
Landau's model the particles do not jump right out into phase space (which leads 
to too many heavy particles in Fermi's picture), but undergo an expansion phase 
before breaking up. The force responsible for the expansion is large in the longitudinal 
direction (the pressure gradient is mainly in the longitudinal direction because of the 
Lorentz contraction) and provides a natural dynamics for the well-known transverse- 
longitudinal asymmetry of secondary momenta." In the light of present-day discussions, 
it is worthwhile continuing the direct quote: "The detailed calculations are made on 
the basis of the classical relativistic hydrodynamics of a perfect fluid, whose energy- 
momentum tensor T'^'^ is: 



where u^{x)[= {uo,u) with uq = y/1 + u ■ u] is the four-velocity field and e, P are 
the scalar densities of energy and pressure. The hydrodynamical equations are simply 
d^T^y = 0. In order to solve these equations one needs in addition an equation of state, 
taken by Landau to be P = e/3 which is characteristic of black body radiation. It is 
not surprising that this is equivalent to the vanishing of the trace of T^^, = e — 3P." 
The expansion phase is taken as scale-invariant {u = x/t) |161j and adiabatic (conserves 
entropy). "The details of the final stage, in which the fluid breaks up into asymptotic 
states of the system, remain quite obscure, as in the parton model." 

The Landau hydrodynamical model was never popular with particle physicists 
because it seemed not relevant to small systems such as p-p coUisons and because 
it predicted a gaussian rapidity distribution, whereas a flat rapidity distribution — the 
rapidity plateau — had been discovered at the CERN ISR |l(i2j . 

It is interesting to note that Bjorken used Landau hydrodynamics in his seminal 




(24) 
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Figure 10. a) Comparison of dN/dy of 7r+ and tt^ in Au+Au central collisions 
at -sJsnn = 200 GeV and Landau's prediction. |163j b) Ratio of data to Landau's 
prediction for rms width {a) of the distribution |163| c) (right) The invariant cross 
sections for 7r+, tt^, , , p and p for pT = 0.4 GeV/c versus j/iab = y — J/^™-* at 
the CERN-ISR [T?)^ . 



paper on A+A collisions [7S] with an initial condition to assure 1 dimensional 
longitudinal 'boost invariant' expansion (resulting a flat rapidity distribution). The 
Id expansion lasts until the rarefaction front which starts at the outer transverse edge 
of the collision zone propagates at the speed of sound to the center and the whole 
fluid (being 'aware' of the edge) starts expanding 3-dimensionally. Recent results 
from Au+Au collisions at RHIC indicate a relatively more Gaussian than flat rapidity 
distribution for identified tt^ |163j (Fig.lTUj). leading to suggestions that the p-p situation 
be reevaluated both with respect to the rapidity plateau and searches for evidence of 
other hydrodynamic phenomena EHl I164j . 

The applicability of the relativistic hydrodynamics of a "perfect fluid" (zero 
viscosity) to both p-p and A+A collisions means that the issue for the QGP will be 
not whether hydrodynamics is applicable, but rather what exactly are the EOS and 
the initial conditions; and what exactly flows — hadrons, constituent quarks or current 
quarks. 



4.1.2. Hard- scattering It is important to distinguish the soft multi-particle physics, 
with limited transverse momentum corresponding to the ~ 1 fm scale of the nucleon 
radius, which is described by constituent- quarks, thermodynamics and hydrodynamics, 
from the high transverse momentum phenomena {px > 2 GeV/c), due to the hard- 
scattering of point-like current-quarks, which correspond to a very short distance scale ~ 
0.1 fm <^ 1 fm |147j . and contribute <^ 1% of the particles produced. These two different 
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pictures of particle production in p-p collisions, the hard and the soft mechanisms, 
indicate that the "elementary" p-p collisions are actually rather complicated processes. 

As noted above (Sec. the observation of particle production at large transverse 
momentum (pr), well above the extrapolation of the low px exponential behavior (see 
Fig.lTTk) proved that the partons of DIS strongly interacted with each other. The subject 
developed before the discovery of QCD, but QCD eventually provided a quantitative 
description of this phenomenon |165j . 




Figure 11. a) (left) Ed?'a{pT)/d^P at mid-rapidity as a function of in p + p 
collisions UMI. b) (right) ^^{GeVf-^ x Ed^a/d^p vs xt = 2pt/^/^ jESl- 



The overall p-p hard-scattering cross section in "leading logarithm" pQCD is the 
sum over parton reactions a + b ^ c + d (e.g. g + q ^ g + q) at parton-parton center- 
of-mass (cm.) energy v^: 

"^'^ - ^ T.U^i)M^2)^^^^^^\cose*) (25) 



dxidx2dcos9* s 2xiX2 

where fa{xi), fb{x2), are parton distribution functions, the differential probabilities for 
partons a and b to carry momentum fractions Xi and X2 of their respective protons (e.g. 
u{x2)), and where 6* is the scattering angle in the parton-parton cm. system. The 
parton-parton cm. energy squared is s = X1X2S, where ^/s is the cm. energy of the 
p + p collision. The parton-parton cm. system moves with rapidity y = 1/2 ln(a;i/a;2) 
in the p + p cm. system. 

Equation I2S1 gives the Pt spectrum of outgoing parton c, which then fragments into 
hadrons, e.g. vr^. The fragmentation function D'f [z] is the probability for a 7r° to carry 
a fraction z = p'"°/p'^ of the momentum of outgoing parton c. Equation |2H1 must be 
summed over all subprocesses leading to a 7r° in the final state. In this formulation. 
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fb{x2) and D^{z) represent the "long-distance phenomena" to be determined 
by experiment; while the characteristic subprocess angular distributions, E"''(cos6'*), 
and the coupling constant, ^^(Q^) = 25 1ii(q2/a2) ' fundamental predictions of 

QCD unnmnzmnH]. 

Equation 123 leads to a general 'xy-scaling' form for the invariant cross section of 
high-pr particle production: 

where xt = 2pT/y/s. The cross section has two factors, a function F{xt) {G{xt)) 
which 'scales', i.e. depends only on the ratio of momenta, and a dimensioned factor, 
1/Pt where n equals 4 in lowest-order (LO) calculations, analogous to the 

form of Rutherford Scattering in QED. The structure and fragmentation functions are 
all in the F{xt) {G{xt)) term. Due to higher-order effects such as the running of the 
coupling constant, as{Q^), the evolution of the structure and fragmentation functions, 
etc, n is not a constant but is a function of xt, a/s. Measured values of n{xT, y/s) in 
p + p collisions are between 5 and 8 |17()j . 

The scaling and power-law behavior of hard scattering are evident from the ^/s 
dependence of the pt dependence of the p + p invariant cross sections. This is shown 
for nonidentified charged hadrons, {h~^ + h~)/2, in Fig. ITTk. At low pt < 1 GeV/c the 
cross sections exhibit a "thermal" exp {—Qpr) dependence, which is largely independent 
of ^/s, while at high px there is a power-law tail, due to hard scattering, which depends 
strongly on ^/s. The characteristic variation with y/s at high pt is produced by the 
fundamental power-law and scaling dependence of Eq. 123 1221 This is best illustrated 
by a plot of 

v^'^^^"'^^ X = G{xt) , (27) 

as a function of xt, with n{xT, ^/s) = 6.3, which is valid for the xt range of the present 
RHIC measurements fFig.lllb). The data show an asymptotic power law with increasing 
Xt- Data at a given y/s fall below the asymptote at successively lower values of xt with 
increasing ^/s, corresponding to the transition region from hard to soft physics in the 
Pt region of about 2 GeV/c. 



5. "A new state of matter" vs. "The perfect fluid". 

With the terminology, concepts and other basics thoroughly discussed, it should be 
straightforward to present the relevant measurements and to evaluate whether the data 
support the conclusions. 

5.1. Particle production, thermal/ chemical equilibrium, strangeness enhancement 

At RHIC energies, as shown in Fig. IHl the (pt) of tt^, K^, p, p, increase smoothly from 
peripheral to central Au-|-Au collisions, and as in pp collisions increase with increasing 
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mass as would be expected for a thermal distribution fEq.lT^. The m-r-scaling property, 
of thermal distributions is well illustrated at SPS energies jl71j where the effects of 
hard-scattering are small (see Fig. IT^). The increase of the (pt) with centrality for all 
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Figure 12. a) (left) Inverse slope of the rriT distribution (colloquially called T) for 
TT*, K^, p, p in 158 A GeV Pb+Pb central collisions |171| . b) (right) Inviarant yield 
of p and p as a function of centrality scaled (by the number of binary collisions iVcoii) 
to lie on top of each other for pt^2 GeV/c | 172| 



particles, nearly linearly with rest mass, is evidence for collective motion ('radial flow', 
see below) and is seen in RHI collisions at AGS [HS], SPS [HI] and RHIC [10111111121113 

energies. The effect is primarily at low px where the slope flattens with increasing 
centrahty as illustrated in Fig. [T2b- 

The semi-inclusive ratios of different particle abundances also vary smoothly as a 
function of centrality in Au+Au collisions at RHIC (see Fig. [T^ ) with a considerably 
larger increase in production than p, p production relative to vr^. However the p/7i~^ 
and p/vr" ratios as a function of pt (Fig. llHb) show a dramatic increase as a function of 
centrality at RHIC |172j which was totally unexpected and is still not fully understood 
(see below). 

The ratios of particle abundances (which are dominated by low p-r particles) for 
central Au+Au collisions at RHIC, even for strange and multi-strange particles, are well 
described (Fig. [H^) by fits to a thermal distribution. 
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with similar expressions for strange particles, where fiB and fis are chemical potentials 
associated with each conserved quantity: baryon number (/x^) and strangeness (fis)- 
This should not be very surprising as the particle abundances in A+A collisions at SPS 
and AGS energies |176j and in p-p jl77j and e~^e~ collisions |178j are also well described 
by the same thermal model "(albeit, only by including a strangeness undersaturation 
factor, 7s < 1) in p-p, e~^e~ and p-|-A collisions, where thermal and chemical equilibrium 
are thought not to be achieved [13)" ? to such an extent that the chemical freezeout 
temperature Tch as a function of fiB (from Eq. |2H1) could be derived, which looks 
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Figure 13. a) (left) K^/tt^ and p^/n^ ratios as a function of ccntrality (A^part) in 
Y^saT/v Au+Au collisions. |63| . b) (right) p/tt and p/n ratio as a function of px from 
the same data |63l 1172) . Open (filled) points are for tt^ (tJ""), respectively. 
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Figure 14. a) (left) Ratios of pT-integrated mid-rapidity yields for different hadron 
species measured in Au-f Au central central collisions at ^/sWn = 200 GeV |^. The 
variation of 7s as a function of ccntrality is shown as an inset including measurements 
for 200 GeV p-p collisions, b) (right) Tch versus /is for thermal model fits as a function 
of ^/snn |174| . where the dashed lines represent the possible phase boundary of the 
QGP transition and the dotted lines represent the result of |175j 



suspiciously like a phase diagram (Fig. Eb)- Of course, as the thermal equilibrium 
properties of the QGP are the subject of interest, it is important to understand how 
or if the thermal properties of the observed hadrons relate to the thermal properties of 
quarks and gluons in the QGP before getting too excited. Also, while these thermal 
models appear to be simple they have many technical details which are beyond the scope 
of this review |179j . 

This brings us to the perplexing question that if all particle ratios are explained by 
a thermal model then how could CERN claim from the strange particle abundances that 
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"all attempts to explain them using established particle interactions have failed" ? This 
requires a somewhat more detailed examination of strange particle production. The 
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Figure 15. a) (left) E802 K/tt vs y in p+A and Si+Au (central) collisions at 
^/snn = 5.4 GeV |180j : b) (right) Integrated yield over all rapidity of per 
participant in Si+Au {y/sNN ~ 5.4 GeV) and Au+Au coUisons {y/sNN ~ 4.7 GeV) at 
AGS [IHIl 



strange particle production at the AGS in p+A and Si+A collisions at 14.6^4 GeV/c 
(yj^ = 5.4 GeV) HHO] and Au+Au collisions at 11. lA GeV/c = 4.7 GeV) [M] 

is shown in Fig. [T3k . In contrast to the pion yield which is roughly constant from 
p+Be to p+Au (not shown), the yield in p+A collisions increases substantially 
with increasing target mass A and and increases in moving from the projectile to target 
region. Another striking feature of the data is that the tt~^/p and K~^/p ratios both 
exhibit a strong rapidity dependence, but the K~^/p ratio exhibits a remarkable target 
mass independence while the 7T~^/p ratios do not. As the authors noted, " The result that 
the / ratio increases with the target mass is somewhat surprising. At AGS energies 
and below, the yield in p+p interactions decreases faster with decreasing energy 
than the vr"*" yield. In a naive picture of successive collisions of the projectile proton 
with target nucleons in heavy targets, a decreasing ratio is expected." In other words, 
strange particles are enhanced in p+A collisions and hadronic models can't explain this. 
In Si+Au collisions, the enhancement of K^/tt^ follows the same trend with rapidity 
but is further enhanced by a factor of ~ 1.5. 

The dependence on centrality and species in Si+Au and Au+Au collisions at the 
AGS is also very striking (Fig. ITM)). The enhancement in the number of per 
participant in Si+A collisions at ^/snn = 5.4 GeV turns on quickly and for Si+Au 
saturates at a moderate centrality. In Au+Au at ^/sNN = 4.7 GeV the number of 
iiT"^ /participant increases slowly and smoothly as a function of centrality reaching an 
enhancement (compared to N-N collisions at the same ^^snn) the same or larger than 
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in central Si+Au, a factor of 3-4 enhancement per participant for compared to N-N 
collisions |181j . The and K~ yields track each other with centrality; and although 
the overall number of compared to K~ increases from Si+Al to Si+Au to Au+Au, 
the K~^/K~ ratio stays constant as a function of centrality for each system |182j . The 
other striking feature of this measurement is that no simple description of the yields 
is given in terms of the number of participants. This tends to go against the simplest 
thermodynamical description, in which the volume (proportional to the number of 
participants) is the only extensive quantity in the problem — this is the thermodynamical 
(rather than quantum mechanical) explanation of why the wounded nucleon model 
works. 

Interestingly, the systematics of kaon production at SPS energies show exactly the 
same effect (Fig. IT^). The total yield of divided by the total yield of 7r+ which is 
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Figure 16. a) (left) Total number of per 7r+ (equivalent to iVpart) in S+S and 
Au+Au coUisons at SPS [1831 15^ : b) (right) dN/dy at midrapidity as a function of 



part 



for strange and multi-strange baryons |184l I185| 



equivalent to /participant, since the pions follow the wounded nucleon model, shows 
the same behavior at SPS energies as at the AGS, which was noted by the authors when 
this data was first presented at the Quark Matter conference in 1999 |183j . However, 
some influential theorists |il86j at this meeting did not seem to be aware of the AGS 
data and this situation persisted in the "New State of Matter" announcement 
The measurements of strange and multi-strange baryons at mid-rapidity were unique 
at the SPS in 1999 since no such measurements had been done at the AGS. They 
showed a progressive increase in enhancement by roughly a factor of 3 per unit of 
strangeness, so that the f2's are enhanced by a factor of roughly 20. Due to the lack 
of measurements in peripheral collisions, the approximate flatness of the enhancement 
for iVpart > 100 |184j attracted some attention as evidence of a pre-hadronic state |186j 
(although the more recent data shown in Fig. ITHb |185j don't look as fiat) and the same 
effect is seen with in Pb+Pb at CERN and Si+Au at the AGS. One pre-hadronic 
state that was neglected was the quantum-mechanical excited nucleon responsible for 
the Wounded Nucleon Model for pions. Multiple excitation of the incident nucleons by 
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successive collisions has little effect on the production of pions, but has a huge effect 
on the production of if it raises the excitation above the K~^A threshold. Also, this 
happens in the initial collisions not in the hadron gas which takes place much later in 
the evolution of the system, so may not be well accounted for in 'hadronic' models. To 
quote Ref. jl86| I187j "The observed strangeness enhancement pattern thus cannot be 
generated by hadronic final state interactions, but must be put in at the beginning of 
the hadronic stage." 

So, does the fact that "models based on hadronic interaction mechanisms have 
consistently failed to simultaneously explain the wealth of accumulated data" jnHl — 
notably the strangeness enhancement both at AGS and SPS — imply that the QGP exists 
at the SPS (and the AGS)? Obviously not! The fact that certain 'hadronic models' do or 
do not explain the data is a statement about the validity of these models (which already 
fail for strangeness in p+A production at the AGS |.18(J|) and is certainly no justification 
for drawing any conclusion about the Quark Gluon Plamsa. The SPS strangeness data, 
either at that time or at present |188j . never supported a QGP conclusion even though 
strangeness enhancement was one of the original proposed signatures for the QGP 
fSec. l3.2rT|) . Strangeness is enhanced in nuclear collisions at all cm. energies |189U19U] . 
even in p+A collisions. Hence strangeness enhancement is an ordinary feature of nuclear 
collisions, not a unique feature of the QGP. 

Interestingly, strangeness production and the region between the SPS and AGS 
energies has drawn renewed interest lately I191j . The integrated yield (over all 
phase space, rapidity, azimuth, pt) is shown in Fig. [T7k fBU] . As noted above, there is 
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Figure 17. a) (left) Ratio of total yields of K^/n^ and /tt^ as a function of 
y/SNN EDI; b) (right) Ratio of ss to uu and dd quarks, A |154j . as a function of -^/s]vjv 
from Ref. |192| with predictions from Ref. |193j . The dotted line is the effect of the T 
variation, with fiB = 0, and the solid line includes the fiB variation. 



an enrichment of in regions of large baryon density due to associated production via 
e.g pp — > pK^A. At the AGS, where the net proton density (/xs) is high at midrapidity 
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(recall Fig. ^ the number of strongly exceeds the number of K~ , which have 
no baryon-associated production channel and are mostly produced in a K~K^ pair. 
The pair production channel increases smoothly with ^/snn while the production 
exhibits a maximum at ^/snn ~ 7 GeV. This maximum is explained |192[ I193j (see 
Fig. ITTb ) by the interplay between the natural rise in strangeness production (e.g. 
K~ , Fig. IT7k ) with increasing ^/sNN (or T) combined with the reduction of associated 
production due to the reduction of fis with increasing T (Fig. IT^ as the nuclear 
transparency reduces with increasing ^/snn and the valence protons move away from 
midrapidity (Fig. El). Of course, the burning issue at this time is whether or not the 
maximum (colloquially the 'horn' [H2]) at ^/sNN ~ 7 GeV is a discontinuity (a real 
signature of a phase transition) or just an example of the very interesting but not 
QGP physics of A+A collisions. Clearly, more precise measurements in this region are 
required. 

5.2. Flow 

A distinguishing feature of A+A collisions compared to either p-p or p+A collisions is 
the collective flow observed. This effect is seen over the full range of energies studied 
in heavy ion collisions, from incident kinetic energy of lOOA MeV to cm. energy of 
^/sNN = 200 GeV |194j . Collective flow, or simply flow, is a collective effect which can 
not be obtained from a superposition of independent N-N collisions. It comes in three 
varieties: directed flow, radial flow and elliptical flow. 




/ I V 

Figure 18. a) (left) Almond shaped overlap zone generated just after an A+A 
collision where the incident nuclei are moving along the ±z axis, and the reaction 
plane, which by definition contains the impact parameter vector (along the x axis) 
Thanks to Masashi Kaneta for the figure |195j . b) (right) View of the collision down 
the z axis: (top) spatial distribution (bottom) momentum distribution after elliptic 
flow (^2) develops 
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Immediately after an A+A collision, the overlap region defined by the nuclear 
geometry is almond shaped (see Fig ITH)) with the shortest axis along the impact 
parameter vector. Due to the reaction plane breaking the symmetry of the problem, 
the semi-inclusive single particle spectrum is modified from Eq. El by an expansion in 
harmonics |196j of the azimuthal angle of the particle with respect to the reaction plane, 
(j) — |197j , where the angle of the reaction plane $r is defined to be along the impact 
parameter vector, the x axis in Fig. 



The expansion parameter vi is called the directed fiow and V2 the elliptical fiow. If 
no collective behavior takes place, i.e. the interaction is merely a superposition of 
independent nucleon-nucleon collisions, then the outgoing momentum distribution of 
the particles would be isotropic in azimuth. However, since the leading participating 
nucleons in the forward region +z (Fig. ITSk ) will interact with many other nucleons in 
the "almond" , they will be pushed away from the rest of the participants, into the +x 
direction, while the —z going participants are pushed towards —x. This is what causes 
the directed fiow, fi, which was discovered at the Bevalac P and is clearly sensitive 
to the Equation of State. For instance if one imagines the almond to be composed of 
billiard balls requiring lots of pressure for a small deformation (hard EOS) a larger vi 
would result than if the almond suddenly melts, perhaps turning into a 'perfect fluid', 
with a much softer EOS |198j . 

The same principles apply to V2, the parameter of cos2(</) — $_r), which (unlike 
Vi) doesn't change sign with rapidity, and hence is non-zero at midrapidity. If thermal 
equilibrium is reached, then the pressure gradient is directed mainly along the direction 
of the impact parameter {x axis in Fig. ITHb) and collective flow develops along this 
direction. If all the particles are approximately at rest in the fluid and thus move with 
the fluid velocity, the transverse momentum distribution will reflect the fluid proflle. 
Hence the anisotropic spatial distribution is carried over to an anisotropic momentum 
distribution through the pressure gradient 

It is important to emphasize that the spatial anisotropy turns into an momentum 
anisotropy only if the outgoing particles or partons interact with each other |197j . Thus 
the momentum anisotropy is proportional to the spatial anisotropy of the almond, 
represented by the eccentricity, e = [Ry — Rl) / [Ry + i?^) ~ {Ry — R^)/ {Ry + Rx) , at the 
time (to) of thermalization. This is due to the fact that the mean number of scatterings 
in the transverse plane is different along the x and y axes |199t l2Um I197j . The mean 
number of scatterings is proportional to the particle density, p = (l/TiRxRy) dn/dy 
(similar to Eq. [T3|) times the interaction cross section (a) times the distance traversed: 
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where = y{x'^), Ry = yiv"^)- Hence one test for hydrodynamic evolution 
to plot V2/e as a function of p = (l/TrRj-Ry) dn/dy (see FiglT^. The data follow this 
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Figure 19. a)(left) V2/£ as a function of particle density p = {l/nRxRy) dNch/dy 
as compiled by Ref. |2U5| for midcentral collisions (~ 13 — 34% upper percentile), 
b) (right) V2/£ as a function of pt for midcentral collisions at RHIC as compiled by 
Ref. [HI 



simple scaling law very nicely. The hydro-limits indicated are for full thermalization 
of the system at the value of e given by the initial nuclear geometry of the almond at 
the time of overlap. If the system doesn't thermalize rapidly, the flow tends to vanish 
because the eccentricity reduces as the system expands [03 1201] . The value of the 
hydro limit is also determined by the EOS through the speed of sound and whether 
the state is hadronic, partonic or undergoes a phase transition |2()2| 1203] . The speed 
of sound, = dP/de (e.g. a simple EOS being e = P/cj.), determines how long it 
takes information about the initial spatial anisotropy to propagate to the whole system, 
~ Ry/cs |2U4j . The hydro-limits shown on Fig. are for a particular choice of an 
EOS with no phase transition |2051 162j . However, recent, more thorough calculations 
indicate that the hydro-limit is reached at RHIC but that V2 at the SPS is considerably 
below it |20^ l2U3j . Because perfect fluid hydrodynamics depends on so few parameters, 
there are many simple scaling tests |2()4j (analogous to EqjnOl and Fig. ITT?|l that can be 
performed to test its validity |194j . 

It is also possible that the increase of V2/e with increasing ^/snn (averaged over 
Pt) is due to the harder pt spectra at larger ^/snn- However the pt dependence of f 2/5 
(Fig. ITM )) clearly increases more rapidly with increasing y^J/vTV- 

Since the eccentricity e is much larger for peripheral than for central collisions, the 
dependence of V2 on centrality has a characteristic shape (Fig. I^Uk ) |2U7j . This was 
one of the first publications from RHIC and showed that w 2 was surprisingly large and 
near the hydro-limits. Another surprise |2U9j (Fig. OUb ) |2U8j was that the V2 followed 
the hydro prediction out to p-r ~ 2 GeV/c and then plateaued at a constant value to 
much higher pt- This was one of the principal arguments for the "perfect fluid" because 
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Figure 20. a) (left) V2 as ftmction of centrality measured as a fraction of the total 
number of charged particles 207^. The boxes represent the expected hydro-limit with 
V2/e = 0.19 (lower edge) and 0.25 (upper edge), b) (right) V2 as a function of px for 
identified particles in minimum bias Au+Au collisions at ^/sWn ~ 200 GeV together 
with a hydro calculation | 208| . 

any modest value of viscosity |21(Jj would cause the V2 to decrease towards zero near 
Pt ~ 1-7 GeV/c (Fig.inii)- 




PT(GeV) p,/n (GeV/c) 



Figure 21. a) (left)u2 as a function of pr in a hydro calculation |21()| for mid-central 
collisions for different values of Ts/tq, the "sound attenuation length" which is zero 
for a "perfect fluid" and increases linearly with the viscosity, b) (left) V2/n vs px/n 
for identified particles, where n is the number of constituent quarks |2(J8| 

As hydrodynamics appears to work in both p-p and A-l-A collisions, and collective 
flow is observed in A-l-A collision over the full range of energies studied, a key question 
is what is flowing at RHIC and is it qualitatively different from the flow observed at 
lower ^/sNN ? One interesting proposal in this regard is that the constituent quarks 
flow |211l| . so that the flow should be proportional to the number of constituent quarks 
Ug, in which case f 2/^^ as a function of pr/ng would represent the constituent quark flow 
as a function of constituent quark px and would be universal. Interestingly, the RHIC 



Recent results in relativistic heavy ion collisions 



35 



data fFig.l21b) seem to support this picture, although the fact that the tt~^ + 7r~ deviate 
most from the universal curve should raise some suspicions as the pion is the only particle 
whose mass is much less than that of its constituent quarks. Another striking hint as to 
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Figure 22. a) V2 vs px for (j> K^K^ and other particles as indieated in y^s^vw = 200 
GeV minimum bias Au+Au collisions j212| : b) same data plotted as V2/n vs pr/n 
where n is the number of constituent quarks, c) (right) V2 of non-photonic e"*" + e~ 
from semi-leptonic heavy flavor decay |21Hj . The solid curve is if both charm and light 
quarks flow, while the dashed curve is if only the light quark flows |214| . 



what is flowing at RHIC is given in Fig. 1^ where the meson fPig.O^) |212j and charm 
particles (Fig. |22fc) detected by their large semi-leptonic decay |2ia EIHl EIZl EH 1219] 
exhibit the same V2 as other particles |22()^ I213j indicating for the that the flow is 
partonic because the hadronic interaction cross section of the (f) meson is much smaller 
than for the other hadrons both in the constituent quark model and as measured in 
photoproduction |22H 1222} I223j . For the charm particles, the V2 of the decay electrons 
follows the V2 of the D{cd,cu) mesons |224[ 1214] . but due to their different masses the 
c and u, d quarks have different momenta for the same velocity required for formation 
by coalescence, so the f 2 of the D mesons is reduced to that of the light quark at lower 
Pt if the c quark itself does not flow. The data (Fig. |22b) favor the flow of the c quark, 
but clearly lots more work remains to be done to improve both the measurement and 
the theory. As no data for either or charm flow exist at the SPS it is difficult to know 
whether the observation of what must be partonic flow for these particles at RHIC is 
qualitatively different from flow at the SPS. 

The last variety of flow is radial flow, which was illustrated in Fig. El Particles or 
partons which travel with the flow velocity, m, acquire increasingly larger kinetic energies 
with increasing mass so that the {pt) increases with increasing mass and centrality which 
is manifested by the flattening of the pt distribution at low pt |197t 1225] , an apparent 
increase in the inverse-slope: 

^ ~ Tapparcnt ^T + Kn (ul) . (31) 

While the elliptic flow is only sensitive to the geometry of the flow profile, the radial flow 
is sensitive to the velocity and hence puts a requirement on hydro models to explain 
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the pt spectra for identified particles as well as the observed V2- While some hydro 
models rise to the challenge, there is no hydro model which also explains the size of 
the interaction volume measured by Bose- Einstein correlations (HBT) |226j . though all 
agree that the flow velocity {u±) ~ 0.5 (in units of the speed of light). Also there is no 
definitive explanation of the decrease of f 2 away from mid-rapidity [SH IH2] which was 
not discussed here. 

5.2.1. Hydrodynamics and Bose-Einstein correlations When two identical bosons 
(usually 7r+7r+ or tt^tt^) occupy nearly the same coordinates in phase space, constructive 
interference occurs due to symmetry of the wave-function imposed by Bose-Einstein 
statistics, a quantum-mechanical effect. This is measured with relatively soft particles 
(Pt < 1 GeV/c) using a two-particle correlation function, C^^ \i27\, which represents 
the probability of detecting two identical pions with momenta pi, Q = Pi ~ P2, 
divided by the product of the probabilities for detecting the individual pions (which is 
the probability of detecting two such pions if they are uncorrelated) . For a classical 
chaotic source, Cf^ = 1, while for a quantum source of identical bosons, Cf^ = 2 |228j . 
For a gaussian source, Cf^ takes the form j228l I229j : 

r>BE( \ _ '^iPliP'i) _ 1 1 \ /' I?2 2 p2 2 p2 2 \(oo\ 

^2 \1) — 'Pjp^ypj^ — -L + ^6Xp|^ — Jtgjdgggidg — -Koutgout ~ -n-longftongJl'J^J 

where the coordinate system is chosen so that giong is along the beam direction, gout 
is parallel to the average transverse momentum of the pair, = [pri +Pt2)/2, and 
gside is the other transverse component, perpendicular to kx- The parameter A measures 
the strength of the correlation and would have the value A = 1 for a fully chaotic 
quantum source. The extracted source dimensions are commonly called HBT radii, 
after a similar technique developed by radio astronomers to measure the angular size of 
radio sources |23()j . 

For a time-dependent source, the energy difference, qo = Ei — E2, of the two pions 
must be taken into account. The radii then become functions of fc^; and -Rout acquires 
an additional dependence on the time duration of particle emission, r, such that for a 
cylindrical source, i?out = -^side + Pr'^'^y where Pt is the average velocity of the pair, 
corresponding to kr- It is important to note that dynamical effects due to final state 
interactions can be large, which makes the interpretation of these measurements a very 
specialized subject [231J. One such example is the case of a source which expands due 
to collective flow. The HBT radii do not increase as a function of velocity but instead 
decrease, since they correspond to "lengths of homogeniety" , or regions of the source 
which emit particles of similar momentum |232j . 

HBT measurements in search of the QGP were primarily motivated by theoretical 
predictions of a large source size and/or a long duration of particle emission which 
would result from the presence of a long-lived mixed phase during a first-order phase 
transition from a QGP to a gas of hadrons |63 | I229[ l233[ I234j . This would be indicated 
by large values of i?out/-Rside > 1.5 — 3. However, beautiful measurements of Bose- 
Einstein correlations with pions at RHIC |235t I239j show no such effect (see Fig.l^ jSSj . 
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Figure 23. The fcr dependence of HBT radii ^ for PHENIX gSH 0-30% most 
central, and STAR 2?9] . as labeled, in Au+Au collisions at y^SNN — 200 GeV 



Although a clear decrease of all the HBT radii with increasing krp is seen, indicating an 
expanding source, -RoutZ-Rside hovers close to or below 1, in disagreement with the hydro 
models [Mlffl- This is the "RHIC HBT puzzle" [22^1 ■ 

Does the HBT puzzle indicate a problem with hydrodynamics calculations, or the 
absence of a first-order phase transition, or both? With regard to the former, hydro- 
inspired parameterized fits seem to be able to explain the data |164j . However, the 
exact solution of 3+1 dimensional relativistic hydrodynamics, using the correct initial 
conditions, EOS, and realistic hadronization is a non-trivial, non-linear mathematical 
problem which has not yet been solved |i5E|- Clearly, the RHIC HBT puzzle is a message 
to the theoretical community that the time has come to make the effort. Regarding the 
absence of a first-order phase transition, this seems to be consistent with the latest 
thinking (see Fig. ^) as well as with measurements of event-by-event fluctuations. 

Measurements of the distribution of the event-by-event average transverse 
momentum of charged particles provide severely small limits on possible non-random 
fluctuations due to a phase transition. For events with n detected charged particles 
with magnitudes of transverse momenta, pTi, the event-by-event average pt, denoted 
Mpj, 1237] . is defined as: 

1 " 

MpT=W=-T.PT, ■ (33) 

i=l 
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A typical Mp^ distribution for a measurement in one centrality class of y/s^^ = 200 
GeV Au+Au collisions in PHENIX |238j is shown in Fig. (data points) compared to 
a mixed-event distribution (histogram) which defines a baseline for random fluctuations. 
The difference between the data and the mixed-event random baseline distribution 
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Figure 24. a) (left) Mp^ for Au+Au collisions at ^/sWn = 200 GeV with 30- 
35% centrality ;238i : measured events (■ ), mixed events (histogram), b) (right) 
Compilation [2111 of Sp^ (in per-cent) for SPS (Pb+Au) and RHIC (Au+Au) 
experiments, as indicated, as a function of ^snn and fiB 



is barely visible to the naked eye. The non-random fluctuation is quantified by the 
difference between the measured and random values of the variance divided by the 
squared mean: 

\ / random 

where /i = (Mp^) and = (^Mp^'^ — {Mp^)"^. If the the entire non-random 

fluctuation were due to fluctuations of the temperature (T) of the initial state |24()j . 
with r.m.s. variation relative to the mean, (Tt/(T), then Ep^, = ax/ (T). FigQ^ |241j 
shows measurements of Ep^ for central collisions at SPS (Pb+Au) and RHIC (Au-f-Au) 
energies, with similar small values (~ 1%) observed in all cases. These results put 
stringent limits on the critical fluctuations that were expected for a sharp phase 
transition, both at SPS energies and at RHIC, but are consistent with the expectation 
from lattice QCD that the transition is a smooth crossover 



5.2.2. Summary and conclusions on flow In summary, collective flow is observed in 
A+A collisions at all cm. energies; but, uniquely at RHIC, hydrodynamics with full 
thermalization appears to describe the V2 and spectra, with constituent quarks as the 
flowing objects. Since hydrodynamics also provides a decent description of p-p data 
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much experimental and theoretical work remains to be done before the smoking gun of 
the QGP can be claimed from the success of hydrodynamics in Au+Au collisions at 
RHIC. Suggestions for further study include measurements of V2 with U+U collisions, 
where the highly deformed Uranium nuclei should provide a strong V2 signal even for 
central collisions |2U3j and the search for collective effects in p-p collisions, where due 
to the difficulty of defining the reaction plane, the search for radial fiow in central (high 
multiplicity) p-p collisions might be the most promising. Also, there is still no unified 
hydrodynamic description at RHIC of V2, pr spectra and spatial sizes determined by 
HBT. Even if it should turn out that constituent quarks fiow at RHIC but not at the 
SPS, one must heed the warning of Ref. [^2]; "it must be kept in mind that constituent 
quarks are not partons: they are effective degrees of freedom normally associated with 
chiral symmetry breaking and confinement, rather than with the deconfinement of a 
QGP." Nevertheless, the description of the medium produced at RHIC as a "perfect 
fiuid" |M1 1242j seems appropriate. 

5.3. Jet Quenching 

One of the major, arguably the major discovery at RHIC, was the observation of jet 
quenching |243^ I244j (Sec. 13. 4j) by the suppression of 7r° and non-identified charged 
hadrons in Au-t-Au collisions at mid-rapidity for large transverse momenta, > 2 
GeV/c. In p-p collisions, particles with J9t > 2 GeV/c at mid-rapidity (perpendicular 
to the collision axis) are produced from states with two roughly back-to-back jets which 
are the result of hard-scattering of the constituents of the nucleon (current-quarks and 
gluons) as described by QCD (Sec. I4.1.2|l . The suppression of high pt particles was 
observed by all four RHIC experiments jHSl IHOl lEH 1^2] and is well cahbrated by using 
measurements in p-p, d-l-Au and Au+Au collisions at ^/sNN = 200 GeV. The pt spectra 
of 7r° from p-p colhsions at i/i = 200 GeV |249j and Au-|-Au collisions at ^/snn = 200 
GeV |250j are shown in Fig. 123 Note that the p-p measurements for Pt > 3 GeV/c are 
in excellent agreement with a pQCD calculation |253j . 

Since hard scattering is point-like, with distance scale 1/pt < 0.1 fm, the cross 
section in p+A (B+A) collisions, compared to p-p, should be simply proportional to 
the relative number of possible point-like encounters |251j , a factor of A (BA) for p+A 
(B+A) minimum bias collisions. For semi-inclusive reactions in centrality class / at 
impact parameter b, the scaling is proportional to TAsib), the overlap integral of the 
nuclear thickness functions |252j . where (Tab) f averaged over the centrality class is: 




and where (Ncou) / is the average number of binary nucleon-nucleon inelastic collisions, 
with cross section aNN, in the centrality class /. This leads to the description of the 
scaling for point-like processes as binary-collision (or Ncoii) scaling. This description 
is convenient, but confusing, because the scaling has nothing to do with the inelastic 





(35) 
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Pt (GeV/c) Pt (GeV/c) 

Figure 25. a) (left) Invariant yield of tt" in Au+Au peripheral collisions at ^/sNN = 
200 GeV (stars) and compared to p-p collisions, scaled by Ncou (circles) |25UI 1249] , 
The black line is a pQCD calculation. |25;-{| The yellow band around the scaled p- 
p points represents the overall normalization uncertainties, b) (right) The same for 
central collisions. 



hadronic collision probability, it is proportional only to the geometrical factor (Tab) f 
(Eq.EH). 

Effects of the nuclear medium, either in the initial or final state, may modify the 
point-like scaling. This is shown rather dramatically in Fig. |231 where for peripheral 
collisions. Fig. l25k . the AuAu data are in excellent agreement with the point-like scaled 
p-p data while for central collisions. Fig. I^Hb . the Au-|-Au data are suppressed relative 
to the scaled p-p data by a factor of ~ 4 — 5 for > 3 GeV/c. A quantitative evaluation 
of the suppression is made using the "nuclear modification factor", Rab, the ratio of 
the measured semi-inclusive yield to the point-like scaled p-p cross section: 

R = ^^^^ = ^^^B 

{TAB)fy<da^j, {N,ou)fy<dN§^ ^ ^ 

where dN^^ is the differential yield of a point-like process P in an A + B collision and 
da^j^ is the cross section of P in an A^A^ (usually p-p) collision. For point-like scaling, 
Rab = 1- 

The suppression, Raa, of 7r° in central Au-|-Au collisions at RHIC (see Fig. EBk). 
although quite dramatic in its own right, is even more dramatic when compared to 
previous data. All previous measurements of nuclear effects at high pt > 2 GeV/c 
in p+A and A+A collisions at lower ^/SNN have given results which are larger than 
point-like scaling (Fig. 126b). a situation called the 'Cronin Effect' |256j and thought to 
be due to the multiple scattering of the incident partons in the nuclear matter before 
the hard-colhsion |257t I258j . The suppression observed at RHIC is a totally new effect. 

Naturally, the first question asked about the RHIC suppression was whether 
it is an initial state effect, produced, for instance, by 'shadowing' of the structure 
functions in nuclei, or a final state effect produced by the medium. Although originally 
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Figure 26. a) (left) Nuclear modification factor Raa{pt) for tt" in central • and 
peripheral ■ Au+Au at y/s^N = 200 GeV. ESDI- The error bars include all point-to- 
point experimental (p-p, Au+Au) errors. The shaded bands represent the fractional 
systematic uncertainties for each centrality which can move all the points at that 
centrality up and down together, b) (right) Compilation of RAAiPr) for production 
in A-l-A collisions from Refs. E 



answered by all 4 RHIC experiments by the observation of no suppression in d-l-Au 
collisions inOl IMl 1121 a clearer answer came from later measurements of QCD 
hard-photon production {gq — > •yq) |259j . Fig. ITfk . and the total yield of charm 
particles {gg cc) deduced from measurements of non- photonic e± (Fig. EZb) 

in Au-|-Au collisions. Both these reactions are sensitive to the same initial state partons 
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Figure 27. a) (left) Raa for pr > 6 GeV/c for prompt photon • and 7r° O production 
as a function of centrality (iVpart) in Au-|-Au colhsions at ^/sNN — 200 GeV |259j . 
Experimental systematic errors for prompt photon are shown by dashes while the 
shaded region on the right indicates the theoretical scale error from the NLO calculation 
of direct photons in p-p |259j . b) Non-photonic electron yield (0.8 < px < 4.0 GeV/c), 
dominated by semi-leptonic charm decays, measured at y^sjvjv = 200 GeV for Au+Au 
(scaled by A^cou) and p-p collisions as a function of centrality (-/Vcoii)- The right-hand 
scale shows the corresponding electron cross section per N — N collision in the above 
Pt range |260| . 



as 7r° production, but the photons, which are themselves elementary constituents which 
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participate in and emerge directly from the hard-scattering, do not interact with the 
final state medium; and the number of cc pairs produced does not depend on whether 
or not the c or c later interact with the medium. Hence the fact that Raa = 1 as a 
function of centrality for these two reactions, which is dramatically different from the 
suppression of 7r°, indicates that the 7r° suppression is produced by the interaction of the 
outgoing hard-scattered parton with the medium, losing energy before it fragments into 
a 7r°. The energy loss is indicated by the shift to lower pj- of the spectrum (Fig. 03b) for 
a given yield. Another possiblity is that the energy loss is hadronic rather than partonic 
(to be discussed below). 



5.3.1. More Surprises It turned out that the d+Au measurement was less clear- 
cut than sketched out in Fig. IHl For particles detected at 90° from the collision 
axis (at mid-rapidity, = 0), the d+Au results showed no suppression, possibly a 
Cronin enhancement |261j (see Fig. OHk). However for more forward emission, larger 
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Figure 28. a) (top) Central • and semicentral O Rcp{pt) ratios for charged hadrons 
at pseudorapidities rj — 0,1.0.2.2,3.2 in d+Au collisions at y^SNN ~ 200 GeV |2()lj . 
b) (bottom) i?(jAu(PT) for same data. Statistical errors are shown with error bars. 
Systematic errors are shown with shaded boxes around the data points. The shaded 
band around unity is the systematic error on iVcoU- 



rj, corresponding to lower values of x (Eq. Rcp, the ratio of binary-scaled semi- 
inclusive yields of central compared to peripheral collisions, showed a huge suppression. 
This was exactly what was predicted if the initial state of the Au nuclei at RHIC were 
a color glass condensate (CGC), in which the gluon structure function at small x is not 
an incoherent superposition of nucleon structure functions but is limited with increasing 
A by non linear gluon-gluon fusion {gg — > g) resulting from the overlap of gluons from 
several nucleons |262j . A CGC can produce either a Cronin-like enhancement or a 
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suppression depending on the initial conditions and the kinematic range covered |263j . 
However, the suppression, if any, at large rj is much less pronounced when a measured 
p-p reference is used and -Rd+Au is plotted fFig. I^Hb). 

There are two lessons from this example, (i) It is important to obtain a p-p 
reference spectrum, since Rcp just shows the relative change of A+A collisions from 
more peripheral to more central which may be confused by a Cronin effect for peripheral 
collisions, (ii) There may be yet another "new state of nuclear matter", the CGC, to 
deal with. This may be particularly important at the LHC where all the physics moves 
to much lower values of x. 

The d+Au data was not the only surprise at RHIC. An even bigger, and totally 
unpredicted, surprise was revealed with the measurement of Raa for identified particles 
other than the 7r° in the range 2 < < 5 GeV/c (Fig. 1^^ . Charm particles, as 
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Figure 29. a) (left) Raa{pt) for non-photonic electrons in central collisions (0-10% 
upper percentile) of Au+Au at ^/snn ~ 200 GeV |213| . b) (right) Rcp for 7r° and 
(j) mesons and ioi p + p in Au+Au at ^/SNN — 200 GeV I264| . The vertical dotted 
bar on the right represents the systematic error in A'coii- The shaded solid bar around 
Rcp — 1 represents the systematic error by which the can move with respect to the 
other data. 



measured by non-photonic + e~ , exhibit a suppression comparable to vr" (Fig. I^Hk). 
indicating a surprisingly strong interaction with the medium |213j . This is is a very 
recent result and has caused quite a commotion since the heavy c quark should lose 
much less energy to the medium than light quarks and gluons j265j . Another totally 
surprising result is that the protons and anti-protons (Fig. EHb) are not suppressed at 
all, but follow point-like scaling. (The Ncoii scaling of the pt spectra can be seen in 
Fig.inb). 

These two results show that particle identification is crucial to unraveling all aspects 
of the physics of heavy ion collisions, even for hard-scattering. In particular, the early 
measurements |243j showed a difference of suppression between tt" and non- identified 
charged particles for 1.5 < Pt < 4 GeV/c (see Fig.lHHk) which led to confusion. However, 
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that effect is now understood to be due to the proton 'anomaly', which appears to 
vanish for pr > Q GeV/c where the Raa from non-identified charged particles comes 
into agreement with the value of Raa for tt", which, notably, remains essentially fiat at 
a value of Raa ~ 0.2 from 3 < < 20 GeV/c. Clearly, something is going on in the 
region 2 < < 6 GeV/c (which has come to be called "intermediate Pt") which needs 
to be understood. 
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Figure 30. a) (left) Raa{pt) for tt^ and non-identified charged hadrons in Au+Au 
central collisions at y/s^N = 200 GeV |266| . b) (right) Raa{pt) for 7r°, prompt 7 and 
77 in the same reaction. |267| 



One interesting proposal to explain the proton anomaly concerns coalescence of 
constituent quarks [2681 12691 l27Uj . If the light quarks have a thermal distribution, 
exp(— then production of protons via coalesence in the intermediate pt region is 
favored over fragmentation from a power-law spectrum. Also, the natural p/ir ratio is 
close to 1 because a proton of a given p^ is produced by the coalescence of 3 quarks of 
Pt/^1 so there is no penalty compared to forming a pion from two quarks of pt/2 because 
[exp(— 6pr/3)]^ = [exp(— 6pt/2)]^. This model nicely predicted that fragmentation 
would become dominant above ~ 5 GeV/c so that the p/vr ratio would return to its 
normal fragmentation value (Fig. Eb) and that the suppression would become the 
same for tt" and p, as nicely shown by the data fFig. IHUk). The model also predicts, as 
observed |264j (Fig. EHb), that the meson (ss) would act like a 7r° meson, which has 
a similar constituent quark configuration, rather than like a proton which has a similar 
mass but a different quark configuration |263]. Similarly, Raa for the rj meson |266j . 
nicely tracks the tt" (Fig. IHUb). The model has an added attraction which is worth 
a direct quote [269 , "Finally, our scenario requires the assumption of a thermalized 
partonic phase characterized by an exponential momentum spectrum. Such a phase 
may be appropriately called a quark-gluon plasma." Unfortunately, other facets of the 
data do not fit this model, as the intermediate px protons and mesons both seem to 
come from jet-like configurations rather then from soft processes (see later discussion). 
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5.3.2. Jet properties from two-particle correlations The study of jet properties via two- 
particle correlations, pioneered at the CERN-ISR |271j . is used at RHIC rather than 
reconstructing jets using all particles within a cone of size Ar = -y(Ar7P~T~(A0)^. 
This is because the large multiplicity in A+A collisions (recall Fig. I2b,c) results in a 
huge energy 7r(Ar)^ x ~ 375 GeV in the standard cone, Ar = 1, for Au+Au 

central collisions. Elliptical flow further complicates the jet measurement because the 
large non-jet background is modulated by cos 20, which is comparable in width to the 
jets studied so far. In Fig. l!?Tk . the conditional probability of finding an associated 
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Figure 31. a) (left) Associated charged tracks with 2 < px < PXt GeV/c per trigger 



charged particle with 6 < < 8 GeV/c for central Au+Au collisions at y/s^N = 200 
GeV as a function of the angle A(j> between the tracks in the range |?7| < 1.4 • compared 
to the data in p-p collisions added to the flow modulated Au+Au background |272| : 
b) (right) Associated charged tracks with 1 < pr < 2.5 GeV/c per trigger charged 
particle with 2.5 < pTt < 4.0 GeV/c after subtraction of flow-modulated background. 
The dashed (solid) curves are the distributions that would result from increasing 
(decreasing) the flow modulation by one unit of the systematic error; the dotted curve 
would result from decreasing by two units |273| . Note that only the jet correlation, 
after background subtraction, is shown. 

charged particle with 2 < pt < PTt per trigger charged particle with 6 < pji < 8 GeV/c 
is shown for p-p and Au+Au central collisions, where the p-p data have been added 
to the large flow-modulated Au+Au background. (Note the offset zero.) For the p-p 
data there two peaks, a same-side peak at A0 = where associated particles from the 
jet cluster around the trigger particle and a peak at A0 = vr radians, from the away 
jet. For Au+Au central collisions, the same side peak is virtually identical to that in 
p-p collisions, while the away peak, if any, is masked by the V2 modulation, and, in 
any case, is much smaller than observed in p-p collisions |272j . The 'vanishing' of the 
away jet is consistent with jet quenching in the medium due to energy loss — the away 
parton loses energy, and perhaps stops, so that there are fewer fragments in a given pt 
range. The fact that the number of associated particles in a cone around the trigger 
particle is the same in Au+Au as in p-p collisions is a strong argument against hadronic 
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absorption as the cause of jet quenching |274j . Since all hadrons would be absorbed 
roughly equally, the associated peak would be suppressed as much as the inclusive 
spectrum in a hadronic scenario, which is clearly not seen. The only escape from this 
conclusion is if the partons or hadrons were so strongly absorbed in the medium that 
only jets emitted from the surface were seen. Of course, since pions at mid-rapidity with 
Pt > 1.4 GeV/c, 7 > 10, are not formed before 14fm, due to the quantum-mechanical 
formation time, they fragment outside the medium, even taking account of the flow 
velocity. Thus hadronic absorption in the medium is not possible for pions. 

The away jet reappears if the pt of the associated charged particles is lowered 
to the range 1 < < 2.5 GeV/c for trigger charged particles with pj-^ in the range 
2.5 < pTt < 4.0 GeV/c (FignTTb) \27'A\ . In the most peripheral collisions, the shape 
of the trigger and away jets looks the same as in p-p collisions fFig. IHTk). However with 
increasing centrality, the away jet becomes much wider (Fig. l!TTb) and possibly develops 
a dip at A0 = vr. Since the outgoing partons travel much faster than Cg in the medium, 
a sonic-boom or mach-cone might develop, as suggested by the dip |276j . There are 
many other ideas to explain the apparent dip, not the least of which is to get a better 
understanding of how exactly to extract the flow effect. 

Study of jet correlations in A+A collisions is much more complicated than the 
same subject in p-p collisions and one can expect a long learning curve. However 
even at this early stage, there is one definitive result from jet correlations |277j . in the 
sense that it is adequate to cast serious doubt on the validity of coalescence models 
to explain the anomalous p/n ratio at intermediate pt- Fig. shows the integrated 
associated particle yields/trigger for p-p, d-l-Au and AuAu collisions from correlations, 
as in Fig. IHTb . where the trigger is identified as either a meson or a baryon in the range 
2.5 < PTt < 4.0 GeV/c and the associated particles, in the range 1.7 < pr < 2.5 GeV/c, 
are not identified. The yield of associated particles/per trigger on the near side, from 
the same jet as the trigger hadron, is the same for meson and baryon triggers as a 
function of centrality, except perhaps in the most central bin; and the same effect is 
seen for the away-side yields. The red-dashed curve indicates the expected trigger-side 
conditional yield if all the anomalous protons in Au-|-Au collisions were produced by 
coalescence. This shows that meson and baryons at intermediate pt are produced by 
hard-processes with the same di-jet structure, and not by soft coalescence. A recent 
paper |278j claimed to demonstrate that two-parton correlations of order of 10% would 
be sufficient to move the red-dashed curve of Fig. |22K into agreement with the data for 
the near side correlations. However the away side correlations, which are the same for 
meson and baryon triggers and show a clear hard-scattering di-jet structure, were not 
addressed. Clearly, the baryon anomaly remains a puzzle. 

5.3.3. Detailed tests of the energy loss formalism At face value, suppression of high 
Pt particles at RHIC appears to support the jet quenching description. Unfortunately, 
as noted in Sec. IH.4^ the properties of the medium can not be derived straightforwardly 
from the data but must be extracted using models |137t 11381 1139j . All the models |262| 
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Figure 32. a) (left) Conditional yield per trigger meson (circles), baryon (squares) 
with 2.5 < pt < 4: GeV/c, for associated charged hadrons with 1.7 < pr < 2.5 GeV/c 
integrated within A(p — ±0.94 radian of the trigger (Near Side) or opposite azimuthal 
angle, for Au+Au (full), d+Au (open ) collisions at y/s^N = 200) GeV [2ZZ|. Shaded 
boxes indicate centrality dependent systematic errors. An overall systematic error 
which moves all the points by 12% is not shown, p-p data are shown for non-identified 
charged hadron triggers, b) (right) log- log plot of central 7r° data from Fig. E5b 



EZniEHDlEHIlEHSEHniEnSI agree that the medium in Au+Au central colhsions exhibits 
a color charge density, expressed as the gluon density in rapidity, of dN^/dy ~ 1100, 
corresponding to an energy density e ^ 30 — 10 GeV/fm^ at an initial thermalization 
time to ^ 0.2 — 0.5 fm, with a time averaged (g) ^ 4 — 14 GeV^/fm, all roughly 100 
times that of cold nuclear matter, or roughly 10 times that of a nucleon, probed in a 
cyhnder around the outgoing parton of radius rz) = l//i = 0.4 fm, where /i = 0.5 GeV 
is the 'screening scale' |137j . which corresponds to y/UrUn for the gluon bremsstrahlung. 

Experimentalists cannot simply accept these numbers pronounced by theorists 
drawing curves through various data points, but must have a way to verify the validity 
of the theory directly from the measurements. Since the theory involves the complicated 
interplay of 3 nuclear effects (Cronin-|-shadowing-|-quenching) |284j the best we can do 
at present is to try to test formulas provided by the theorists to approximate their 
results. For instance, although the total energy loss AE in a static medium is supposed 
to be independent of the energy of the parton, E, and proportional to L^ |134j . once 
expansion is added, "the radiative spectrum and the mean energy loss can be related to 
the soft gluon rapidity density" |279j , and {AE) becomes proportional to L but remains 
independent of E apart from logarithms! 



(AE)^ — Lin— - + ■■■ , 37 

where Cr = 4/3 (3) is the Casimir invariant for quarks (gluons) and A± ^ irRxRy is 
the effective area of the collision. 
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The average energy loss (AE) can be measured from the shift in the px spectrum 
at constant Yield |136j (see Figl32b). It is apparent from this log-log plot of the binary- 
scaled p-p spectrum and Au+Au central collision spectrum for vr", that the p-p spectrum 
is a pure power law for > 3 GeV/c, Ed?a /dp^ = Ap^, with n = 8.1 ± 0.1 |285j . and 
that the Au-|-Au spectrum is shifted down in pt from the scaled p-p spectrum by a 
constant fraction, —Apt/pt = >S'ioss, since the spectra are parallel. This explains why 
Raa is a constant (see FiglHUb). since the relation between ^loss and Raa is particularly 
simple for a power law: 

^Ioss(Pt) = 1 - i?AA(pT)'/^"-'^ . (38) 

The dependence of Eq. IHZI on the gluon density can be tested by measurement 
of S'loss as a function of centrality (recall Fig. ITfj) [UHj . The L dependence of Eq. EH 
for constant initial condition can be studied by measurement of S'loss as a function of 
angle with respect to the reaction plane, where the average L of the system can be 
calculated |285j (Fig. l33k). The compilation of S'loss vs L by this method (Fig. l33b) 
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Figure 33. a) (left) S'loss as a function of angle = — with respect to the 
reaction plane in y^SNN = 200GeV Au+Au collisions as a function of centrality ' 285j . 
b) (right) S'loss versus L^, the length from the center of the almond, combining the 
various data points from (a), with different centrality bins represented by different 
symbols. 

shows a reasonably linear behavior with L with the totally striking result that S'loss 
appears to go to zero at finite L ~ 2 fm, suggestive of a formation time. This may 
explain why attempts to explain the azimuthal anisotropy f 2 at large px solely in terms 
of energy loss failed |286| 1287] . Again, it should be noted that this type analysis is in 
its infancy and a long learning curve is to be expected. 



5.3.4- The smoking gun? The jet suppression observed at RHIC is unique in that 
it has never been seen in either p+A collisions or in A-fA collisions at lower ^/s^n 
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and it probes the color charge density of the medium. Many questions and unsolved 
problems remain which are under active investigation, but this effect comes closest of 
all, in the author's opinion, to meeting the criteria for declaring the medium a Quark 
Gluon Plasma: 

• There is no such effect in p+A collisions at any ^/sNN 

• It is not the 'ordinary physics' of A+A collisons since it only occurs for ^ysNN > 30 
GeV 

• In all discussions of the effect, the operative 'charge' is color and the operative 
degrees of freedom are quarks and gluons. 

5.4- J/"^ Suppression 



I 

OJ 

< 

DC 



1.2 
1.1 

1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 



E772, p + A-> 

Integrated Cross Section Ratios 




Ca Fe 



.96 



, .92 



□ Y 



10 100 

Mass Number 



< 



4- 



3 - 



0.5 - 
0.4- 



^ p(450GeV/c)-p,d(NA5i) 
A p(450'GeV/c)-A (A=C,Al,Cu,W) (NA3S) 
* p(200 GeV/c)-A (A=Cu,W,U) (NA38) 
T "'0( 16x2011 GeV/c)-Cii.L) (NA38) 
■ '^S(32x200 GeV/c)-U (NA38) 
""pb(208xl58'GeV/c)-Pb (NA50) 




200 GcV/c 



10 



10 



10 



10 



10 



10 



^projectile ^ tai'get 



Figure 34. a) (left) A dependence of total cross section for ^P, , T production in 800 
GeV p+A collisions |288| . The per nucleon ratio between heavy nuclei and deuterium 
is shown versus the mass number, A. b) (right) Total cross section for J/'^ production 
divided by AB in A+B collisions at 158-200A GeV 



Although J/\l/ suppression was thought to be the 'gold-plated' signature of 
deconfinement since 1986 jl2Uj . there were problems from the very beginning. The 
principal problem is that the J/\l/ is suppressed in p+A collisions (see Fig. l3l|l. Point- 
like processes such as hard-scattering scale relative to p-p collisions like A{AB) for 
p-|-A (A+B) minimum bias collisions, but unlike Drell-Yan which exhibits an A^'^^ 
dependence, the J/\l/ is suppressed by A^-^^ in p+A collisions |28Hj (Fig. l!Mk). The 
suppression continues as (AB)^'^"^ for minimum bias A+B collisions |126j (Fig. 1Mb) at 
the SPS until the heaviest system, Pb+Pb, where the suppression increases by ~ 25%. 
However the suppression is much more impressive as function of centrality (Fig. IH^k). 
and it was claimed that the expected discontinuity |126j due to the phase transition, 
and the sequential suppression of the Xc ^' and direct J/\l/ had been observed |127j . 
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Figure 35. a) (left) Cross section for J/'i divided by Drell-Yan as a function of 
centrality measured by Et for 158 A GeV Pb+Pb collisions PH5) . The solid line 
represents the expected suppression for cold nuclear matter. The plot may be converted 
to Raa by dividing by 34. b) Raa for J/^ production as a function of centrality 
(-^part) at y/sNN = 200 GeV for the systems and rapidity ranges indicated |296| . The 
solid and dashed lines represent the expected suppression in cold nuclear matter for 

Au+Au i2nn! 



To say the least, the interpretation of the suppression in terms of a QGP was 
controversial, and the effect at the SPS could be explained by hadronic |29H I292j and 
even thermal models |298j . The clincher would be the measurement of the J/'^ at 
RHIC. All the initially produced J/^, the ones not suppressed at CERN, would be 
totally suppressed in the much hotter denser QGP at RHIC, which would prove that 
the J/\E' suppression at CERN was indeed the result of deconfinement. However, in the 
ensuing years a new problem developed when it was realized that if a QGP were indeed 
produced, the thermal c, c quarks would recombine |294^ 1295] to form J/\E'. Thus, if the 
J/\l' suppression were the same at RHIC as at CERN, this would imply that RHIC, 
not CERN, had discovered the QGP since all the initially produced J /"if, which were 
suppressed at CERN by whatever mechanism, would be totally suppressed at RHIC, 
leaving only the thermal J/^ produced by recombination in a QGP. The new problem is 
that nobody would believe this explanation. Incredibly, this is exactly what happened 
(see Fig. the J/\E' suppression, expressed as Raa turned out to be the same at 
RHIC [296] as in the famous SPS measurement [289j 

This effect is illustrated quantitatively in Fig. EHl Models EHS EHZ] which 
reproduce the SPS J/\Ef suppression with or without a QGP predict a near total absence 
of J/\E' at RHIC beyond 150 participants without recombination. With recombination 
turned on, the RHIC data are reproduced (Fig. IHUb). with one notable exception |298j 
which predicts a larger Raa for J/\E' than at CERN. In fact, a J/\l/ enhancement would 
have been the smoking gun for the QGP. It will be interesting to see whether this occurs 
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Figure 36. Raa data from Fig. 135b together with predictions from models (a) without, 
(b) with recombination of cc quarks. 

at the LHC. 

Does the agreement of the CERN and RHIC data for Raa in J/"^ production 
ehminate J/\l/ suppression as a signature of deconfinement? Is it possible that the J/\l/ 
(c, c) is no different in its QGP sensitivity than the 0-meson (s, s)? The recent increases 
in the predicted dissociation temperatures |123j . give a possible way out jl25j . Satz has 
proposed that the Xc and the were suppressed both at CERN and at RHIC, but in 
neither place was the direct J/\l/ suppressed, since a temperature ~ 2Tc was not reached. 
However temperature sufficient to melt the J/\l/ should be reached at the LHC. Does 
this mean that we have to wait until 2009 to prove or disprove an idea proposed in 1986? 
Fortunately there are other tests to be made on the RHIC data. If the Xc or were to 
be observed at RHIC, that would settle the issue. For recombination to be true, J/\l/ 
flow should be observed \2W7\ . and both the rapidity and the pt distributions of J/\l/ 
should be much narrower due to recombination than for directly produced J/\l/ j299j . 

One thing is perfectly clear from this discussion: the claim of the QGP discovery 
from J/\l/ suppression at the SPS was, at best, premature. 

5.5. Thermal photons, leptons 

The hardest experiments in relativistic heavy ion physics involve production of lepton 
pairs with low mass and photons at low px due to the huge backgrounds. A measurement 
of e+e~ pairs in 158A GeV/c Pb+Au (Fig. EZK) had created great excitement because 
it looked as if the "extraordinary signal" of a thermal p° dilepton peak |118j had 
appeared, which stimulated much model-making. A new high resolution measurement, 
this year, in 158A GeV In+In collisions, with beautiful u fi'^fi~ and 
peaks |H()2j should help clarify the issue (Fig. ITTb). 

It is generally agreed that the observation of thermal photon or lepton pair 
production would represent the 'black body radiation' of the QGP. An exciting step in 
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Figure 37. a) Invariant mass spectrum of e^e~ pairs in central 158 A GeV/c Pb- 
Au collisions, where each or e~ has px > 200 MeV/c, normalized to the number 
of charged particles in the detector ^00 , compared to a coctail of known hadronic 
sources. The 'enhancement' from 0.2 < < 0.6 GeV occurs for the pair transverse 
momentum p|f < 0.5 GeV/c [80 1| . b) (right) Invariant mass spectrum of from 
minimum bias 158 A GeV/c In+In collisions 302 . Total data (upper histogram), 
combinatorial background (dashed), fake matches (dot-dashed), and net spectrum after 
subtraction. 

this direction was also made this year. Prompt photon production with > 1.5 GeV/c 
had been observed in 158 A GeV Pb+Pb collisions |303j (Fig.lHHk) but these photons were 
consistent with being from hard-scattering, with an enhancement compared to point-like 
scaled p-p collisions. Measurements of prompt photon production by a new method, 
internal conversion, were presented this year from Au+Au collisions at a/sjvtv = 200 
Gev |304j ■ which showed a clear exponential spectrum in the range 1 < pr < 3 GeV/c 
(Fig. l38b ) . Is this the smoking gun of the QGP or do p-p collisions also show a low pt 
exponential for direct photons as they do for pions? Incredibly, this is not known from 
previous measurements and is an example of how unknown aspects of p-p physics are 
unearthed by A-l-A measurements. The p-p measurement is now underway at RHIC. 
This is an example of what may be faced with any apparent A-l-A discovery at the LHC, 
where the p-p physics is totally unknown |3()5j . 

6. Conclusion 

Was a new state of matter discovered in Pb-|-Pb collisions at ^/sNN ~ 17 GeV? Possibly, 
but it surely wasn't the gaseous QGP ^21 IM IM EHU 1123 ES] envisioned at that 
time; although the logic — since the data couldn't be described by hadron gas models, it 
must be something else, and since the only other model (at the time) was the gaseous 
QGP, that's what we must have discovered — deserves, at least, a smile for its audacity. 
Nevertheless, this region of A-l-A collisions is interesting as the transition from "Baryonic 
to mesonic freezeout" |193j . so it is very possible that a new state of matter does exist in 
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Figure 38. a) (left) Invariant yield of prompt photons from 158A GeV Pb+Pb 
collisions together with point-like scaled p+A measurements. Data points with 
downward arrows indicate unbounded 90% Confidence Level upper limits, b) Invariant 
yield of prompt photons as a function of px from central ^/sJtn Au+Au collisions j,''i()4| 
via photon detection O ^ and from measurements • of e'^e" pairs with 90 < < 300 
MeV, where only internal conversions from prompt photons (and some rj — > ^e~^e~) 
contribute, compared to < niee < 30 MeV, where all internal conversions (mostly 
TT*^ — > je^e~) contribute. 



this region, such as an initial state of baryon resonances, raised over the threshold 
by successive excitation, or even something more interesting [H^. Clearly, this region 
deserves further study. 

Was the gaseous QGP discovered at RHIC? Certainly not! The medium discovered 
at RHIC has the properties of a hot, dense, strongly interacting, (possibly) perfect fluid, 
the sQGP \<r)4:\ I242j . Although the phenomena so far observed at RHIC are far from 
being understood, the scope of experiments has moved from the 'discovery' phase to 
the 'exploration' phase, to characterize in detail the properties of the medium. Has 
extended deconfinement [22] over a large volume been discovered at RHIC? There is 
no evidence for (or against) extended deconfinement. However the observation of ~ 10 
times the color charge density of a nucleon in a volume roughly the size of a nucleon, 
where all the color charge is active on a partonic test probe as it passes through the 
medium, is indicative of deconfinement on a small scale. This is not unreasonable, 
since there is nothing yet in the theory of QCD that will reduce the confinement radius 
with increasing temperature and density. Such a collection of quarks and gluons is not 
describable in terms of ordinary color-neutral hadrons, because there is no known self- 
consistent theory of matter composed of ordinary hadrons at the measured densities. 
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The PHENIX experiment jHS] calls this state, "dense partonic matter." The situation 
is reminiscent of the QGP as described by Satz in a previous review article in this 
journal jlH] (Fig. EHl), except that the distance scale in Fig. EHb is comparable to the 
size of a nucleon rather than the size of the whole nucleus. 




1 .6fm 



Figure 39. "Strongly interacting matter as nuclear matter at a density of closely 
packed nucleons (a) and at a much higher density (b)." 021 The dimension of (b) is 
based on RHIC results. 



Is the gaseous QGP ruled out? No. It is possible that, at the LHC, the much larger 
y/SNN = 5500 GeV will produce a medium with an initial temperature high enough so 
that as{T) will become small enough so that the interaction becomes weak enough that 
a gaseous QGP is produced. This is also the region, where another possible new state of 
nuclear matter, the Color Glass Condensate |262j . should become dominant, if it exists. 

Clearly, the phase diagram of hot, compressed nuclear matter is much richer and 
more complex than originally envisioned, and much work remains to be done at all 
temperature and density scales. 
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